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I . INTRODUCTION 


This is an Interim Report on a research study being conducted 
for NASA/Marshall Spaceflight Center by Physical Sciences Inc. on the 
subject of laser-heated thrusters. The objective of this study is the 
development of a computer program useful for the design of the thrust 
chamber for a CW laser-heated thruster. 

Physical Sciences has been studying the physics and fluid mech- 
anics of laser-heated CW thrusters under NASA sponsorship for about six 
years. The first study ^ as reported in Ref. 1, xronsidered methods 
of absorbing laser energy into gases, a simple flxiid-mechanical model 
of a laser-heated gas flowing in a nozzle, and the stability of the flow. 
The second study. Ref. 2, concentrated on a CO^ 10.6 ym laser heating 
a pure hydrogen propellant gas. A detailed model was constructed for 
laser absorption, the flow of dissociating and ionizing hydrogen, and 
the radiation from the hot hydrogen. Methods of initiating the absorp- 
tion of laser energy in time were also investigated, and some study was 
made of the use of particles to cibsorb radiation from the hot hydrogen. 

An important element of a pure hydrogen CW laser-heated thruster 
is the method of initiating the absorption of laser energy in space. 

For pure hydrogen, the only absorption mechanism is inverse Bremsstrahlung, 
which requires the presence of electrons. But electrons are only pro- 
duced in hydrogen by ionization which occurs at temperatures above 7000 
to 11000 K, depending on the pressure. Thus, only after the hydrogen gets 
hot can the laser energy be absorbed. Some other mechanism must operate 
to heat the gas from its inlet temperature to these high temperatures 
where electrons are produced. 

Such a mechanism has been observed to exist in air and a few 
other gases. It has become known as a laser-supported combustion (LSC) 
wave. Such a wave is a region of rapidly increasing temperature in a 
laser-heated gas. The laser energy is absorbed in the hot gas near the 
rear of the wave. This hot gas radiates and conducts energy forward 
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down the temperature gradient to the cooler part of the wave and the 
incoming cool gas. This radiation and conduction then provides the 
mechanism for heating the cold gas up to the temperature where electrons 
are produced and inverse Bremsstrahlung can act to absorb energy and 
heat the gas. Such an LSC wave is the laser-heated analogy of an ordinary 
deflagration wave heated by chemical combustion. It is a low-speed, 
nearly constant pressure wave in which the temperature increases, and 
the density decreases. In order to conserve mass, the velocity increases 
to keep the mass flow per unit area constant. 

This LSC wave provides a possible mechanism for initiating the 
absorption of laser energy in space in a pure hydrogen thruster. The 
wave is stationary at the beginning of the heating region, and as the 
hydrogen flows through it, it heats ^and absorbs the laser energy. 

Although no hydrogen LSC waves have been studied experimentally, they 
have been observed in air (Ref. 3).* They are found to have a charac- 
teristic propagation speed, which depends on the pressure level and laser 
beam intensity. If they are to remain stationary in a flowing gas, the 
incoming gas must flow at this LSC wave speed. 

In the hydrogen 10.6 \m laser-heated thruster concept studied 
in Ref. 2 an LSC wave was postulated as the mechanism for initiating 
the absorption in space. However, no detailed study of the wave was 
performed. A wave speed for a given intensity was assumed based on a 
simple modification of a theory previously presented for air. It was“ 
important, however, to examine LSC waves in hydrogen to determine their 
speed vs. intensity relation- Therefore, in Ref. 4, a detailed model of 
hydrogen LSC waves was constructed, and their speed-intensity character- 
istics were calculated as a function of pressure. The subsequent accel- 
eration of the hot hydrogen emerging from the wave was also considered for 
one particular case. 

The studies described in Refs. 2 and 4 show the main features of 
laser-heated thrusters using pure hydrogen. Because the absorption mech- 
anism is a high temperature one, these thrusters operate at peak temperatures 

* Since the present work was performed, LSC waves have been observed in 
hydrogen by Conrad at the Army High Energy Laser Laboratory 



of about 20000 K. Even after acceleration to a throat, the gas is nearly 
16000 K. These high temperatures lead to the possibility of very high 
specific impulses, perhaps near 4000 s, but they also lead to very large 
convective and radiative heat fluxes to the containing walls, and therefore, 
large energy transfer out of the gas. The case calculated to the throat 
in Ref. 4 had only 58% of the laser energy still in the gas when it 
reached the throat. It was also found that the LSC wave in hydrogen 
had a threshold intensity, below which a wave could not exist. This 
is in agreement with results for air LSC waves. If losses due to radiation 

and sideways convection are included, the threshold intensity increases. 

2 

Threshold intensities are quite high, above 3 E4 W/cm at 10 atm, so 
that for laser power in the kilowatt range the size of the laser spot is 
quite small, leading to a small thruster cross-section. For example, a 
10 kW laser and a pressure of 10 atm would require a spot no larger than 
0.65 cm in diameter. 

In view of the high temperature and small size associated with 
pure hydrogen rockets at moderate laser powers, a study of other alter- 
natives seems warranted, to determine how thrusters can operate with lower 
temperature and specific impulse. One idea is to seed the hydrogen with an 
easily ionizable gas, such as an alkali metal vapor. This seed will yield 
electrons at much lower temperatures than pure hydrogen, and so permit 
inverse Bremsstrahlung absorption to start at lower temperature. This 
lowers the whole temperature level of the thruster, and reduces radia- 
tion and convection out of the gas. An LSC wave mechanism will still be 
needed, however, to heat the gas from the inlet temperature up to the 
temperature at which sufficient ionization has occurred to permit absorp- 
tion of laser energy. But this will be an LSC wave in seeded hydrogen, 
and will have a different wave speed vs. intensity characteristic than 
in pure hydrogen. Such an LSC wave is being studied under the present 
program. 

A second alternative to the pure hydrogen thruster is to seek 
an absorber which operates at the inlet gas temperature, and continues 
absorbing until the alkali metal seed produces electrons to continue 
absorption. Such a low temperature absorber operates by molecular absorption. 
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and must be chosen so as to be "tuned" to the particular laser radiation 
being used. Based on the absorption studies of Ref. 1, a likely can- 
didate for 10.6 )Jm radiation is water vapor, which will absorb this radi- 
ation until it dissociates at around 4000 K. By that temperature, the 
alkali metal seed will have produced enough electrons to let inverse 
Bremsstrahlung take over as the absorber. With both a low temperature 
absorber, H^O, and a medium temperature absorber, alkali metal vapor 
(and a high temperature absorber, hydrogen, if needed) there is no need 
'for the use of the LSC wave mechanism to heat the gas. It is possible 
for the absorption to proceed continuously from the inlet temperature 
to the temperature at which all the laser energy in the beam is absorbed. 
We will refer to this as the continuous heating case, in contrast to the 
LSC wave case, where only medium and high temperature absorbers are used. 

r 

If one is interested in radiation from CO lasers, around 5.3 ym, 
then Ref. 1 suggests that either NO or H^O can function as the low 
temperature absorber. In thxs case, there is also a molecular absorber, 

CO , which survives to medium temperature (at least 6000 K) , and could 
be used instead of the alkali metal vapor seed. It is advantageous not 
to use the latter seed, if possible, because its large molecular weight 
reduces the specific impulse. 

In any of these thruster concepts it is very important to 
consider the radiative and convective heat flux to the walls, since it 
may prove difficult to cool the walls if these fluxes are too high. The 
objective of high specific impulse requires heating the gas to rather 
high temperatures (4000 to 5000 K, or higher) and the resulting heat 
fluxes will be substantial. One of the purposes of the present work 
IS to develop a computer model which will calculate these fluxes, so 
the difficulty of handling them may be evaluated. 

The objective of the present work, as stated above, is to develop 
a computer program useful for the design of a thrust chamber for a CW 
laser-heated thruster. The discussion above indicates that there are 
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two basic modes of gas heating. One is the LSC wave, seeded hydrogen 
mode* The other is the molecular absorber, seeded hydrogen, continuous 
heating mode. These two modes will require two different computer pro- 
grams. The first will be an extension of the pure hydrogen LSC wave 
code developed in Ref. 4, and the other will be an axisymmetric channel 
flow code including laser energy absorption, and convective and radiative 
heating. 

The particular laser wavelength/propellant gas mixture combina- 
tions we have considered are: 

1- CO (10.6 \m) radiation/H - Cs 

2. CO^ (10.6 yro) radiation/ - H^O - Cs 

3. CO (5.3 ym) radiation/ H^-NO-CO or H^-H^O-CO. 

The first is the LSC wave heating mode and the second and third are in the 
continuous heating mode. 

To develop the computer model a number of physical properties 
of the gas mixture must be modeled : 

1. Laser absorption properties. 

2. Radiative and absorptive properties. 

3. Thermodynamic and transport properties. 

These properties must be incorporated in a flow model which includes 
the relevant effects of heating by laser absorption, gas radiation emis- 
sion and absorption, as well as the usual fluid flow effects of convection 
and heat transfer. 

In this interim report we will describe the progress we have 
made in each of these areas. Section II deals with absorption of 10.6 
and 5.3 ym radiation by inverse Bremsstrahlung, and by the molecules 
H^O, CO, and NO. Section III develops a model for the absorption coef- 
ficients of the important radiative species contained in the gas mixtures 
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listed above, from which their radiative and absorptive properties can 
be obtained. Section IV develops a model for the thermodynamic and trans- 
port properties of a mixture of hydrogen and cesixm up to temperatures 
of 20,000 K, for use in the LSC wave heating mode. For the other mixtures, 
in the continuous heating mode, standard chemical equilibrium codes can 
be used up to the temperatures of interest for continuous heating, which 
are approximately 6000 K. In Section V, the continuous heating flow model 
is described. As yet, the extension of the hydrogen LSC wave code to 
hydrogen- cesium mixtures has not been completed, so its description is 
reserved for the final report at the completion of the present work. 

The final section of this report. Section VI, sets forth the conclusions 
whxch can be drawn so far from the work completed. 

Before proceeding to the details of the modeling, it is of in- 
terest to consider the over-all energy, thrust, specific impulse and 
mass flow relations for a laser-heated thruster, or indeed any thruster. 
This will give us some idea of the parameters needed for a desired level 
of performance. There are simple relations between thrust Th, mass flow 
rate m, vacuum specific impulse power in the gas converted to 

exhaust velocity, P^. The exhaust velocity is gl^ , so that 


Th = 

* ■ • 

Plots of m and Th vs I are shown in Fig. 1-1, with an auxiliary scale 
2 

of P /m = (gl ) /2. The dimensions of P /m are energy per unit mass, 

G sp G 

the same as enthalpy, so for a given gas at a given pressure, a temperature 
scale can be drawn on the axis. The parameter on the curves is 

the power in the gas. 

We see, for example, that at I = 2000 s, for a power in the 

sp 

gas of 10 MW, the thrust is 1 kN, and the mass flow about 0.05 kg/s. For 
10 kW in the gas, the thrust is 1 N and the mass flow is 5 E-5 kg/s, or 
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0.05 g/s. The enthalpy P_/in of this gas is 2 e 8 J/kg, which corresponds 
to a temperature of 5000 K in hydrogen at a pressure of 30 atm. 


This simplistic calculation does not account for energy lost 

from the gas by convection and radiation; the laser power must certainly 

be larger than P , and the enthalpy in the heating region must be larger 
G 

than P /m. However, it does give us an idea of how hot the gas must 
G 

be to attain a desired It would appear that reaching 6000 K in 

hydrogen will put enough energy into the gas to reach values of in 
the interesting range of 1000 to 2000 s. To reach 4000 s requires tem- 
peratures of the order of 16000 K. 
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II. LASER ABSORPTION COEFFICIENTS 


In a laser-heated thruster, the gas is heated by energy which it 
absorbs from the incoming laser beam. This energy transfer must be 
accomplished by some mechanism which permits the coupling of the mono- 
chromatic energy in the laser beam into the gas. There are a number of 
such mechanisms. A common one is inverse Bremsstrahlung, which involves 
the interaction of electrons with either ions or neutral particles to 
absorb the laser energy- This mechanism can work for any type of particle 
and any laser beam wavelength, since it does not require a resonance 
between the laser and the gas particles * 

A second possible absorption mechanism is molecular absorption, 
in which the laser energy is absorbed by means of a resonance with some 
transition in the internal state of a particular molecule. This does 
require matching the absorber to the laser energy wavelength. 

Since the inverse Bremsstrahlung mechanism requires electrons, it 
can only operate at temperatures where they are present. In a propellant 
gas like hydrogen, electrons are made only at high temperatures , since the 
hydrogen must first dissociate before it ionizes. Therefore, inverse 
Bremsstrahlung will become an effective absorber only at around 10,000K 
or higher. One can reduce the temperature where electrons appear by in- 
troducing small amounts of an easily ionized gas, like the alkali metals, 
which can ionize at a few thousand degrees. Such a seed gas will permit 
absorption by inverse Bremsstrahliing at much lower temperatures than pure 
hydrogen. But even with a seed, absorption cannot occur at low tempera- 
tures if electrons are needed. If one starts with cold gas (T < 1000 K) , 
some other heating mechanism will be required to heat the gas up to 
temperatures where electrons are made so that inverse Bremsstrahlung 
can take over. 
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One possible low temperatxire heating mechanism is conduction and 
radiative heating. This is the mechanism used in laser-supported 
combustion (LSC) waves, and its application to pure hydrogen is described 
in Ref. 4. It leads to temperatures on the order of 20,000K in the 
hottest region of the gas, and makes the thruster difficult to design 
because of the necessity of handling such hot gases with their attendant 

N 

high rate of loss by radiation and convection. In the present work, 
the effect of seeding with cesium on the hydrogen LSC wave is to be 
explored to see how that moderates the temperature level at which these 
waves operate. 

In addition to the LSC wave mechanism of radiation and conduction, 
another possible low temperature mechanism is molecular absorption. This 
can operate from room temperature up to the temperature at which the 
molecules dissociate, and if the latter is high enough so that seed electrons 
are present, the high temperature inverse Bremsstrahlung mechanism 
can take over. This gives the possibility of continuously heating 
the gas by absorption from the inlet temperature until the laser energy is 
completely absorbed, without the need for an LSC wave to form. This 
continuous heating mode is also to be explored in the present work . 

The laser wavelengths of interest here are lO.Gym and 5.3ym. The 
molecular absorbers considered have been based on the study described 
in Ref. 1. For lO.Gym we have considered and for 5.3ym, NO, 

and CO. 


The absorption of a gas to a laser beam of intensity I is described 

by an absorption coefficient k defined by the equation 

L 

dl/dx = “k I (2-1) 

L 
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describing the rate the beam intensity decays along the direction of' 
propagation. In this section we shall describe the absorption coefficients 
for the various absorbers and absorption mechanisms considered. 

2.1 Inverse Bremsstrahlung Absorption (e + X + hV e + X) 

The inverse Bremsstrahlting mechanism involves interactions between 
electrons, and either ions or neutrals. For ions, the interaction is by 
Coulomb forces, and is independent of the type of ion. For neutrals, it 
depends on the particular neutral particle involved. 


Electron-Ion Inverse Bremsstrahlung (e + I +hv->e+l ) 

Ref. 1 gives the electron-ion absorption coefficient in terms of 
electron and ion number densities n , n as 

ili X 

’"lei = "e 

where A is the laser wavelength, and is the electron-ion absorption 
cross-section, given in Ref. (2-2) as 

4 / 27t V* 

“ei ° 3 \ 3m^kT ] ^^^4 ^ 

Here Ze is the ion charge, e the electron charge, c the speed of light, 
hp Planck’s constant, and the electron mass. For singly charged ions, 
the expression for a__ becomes 

til , 
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=1.37 E-23 


(2.3a) 


'El 


3 3 

^ / T (cm in cgs) 


3 i 5 

= 1.37 E-27 A / T (m in mks) 


(2.3b) 


The expression given so far is a semi-classical one. To account for 
quantum-mechanical effects, Eq. (2.2) is multiplied by a Gaunt factor G, 


which depends on temperature . We have 
where we fitted a parabola to Fig. 5. 
figure areY^ = 157,900/T/ u = b^c/AkT 
ing table from the Figure: 

Y ^ T 

3 52633 

10 15790 

10^ 1579 

We fitted a quadratic to this table to 

G = 1.04 + 3.74E- 


gotten this factor from Ref. 5, 

For A =: lo.oym the parameters of the 
= 1357/T. Then we read the follow- 


u 

9ff = G 

.0258 

2.1 

.0859 

1.55 

.859 

1.1 

get a Gaunt 

factor 

5T - 3.28E 

2 

-lOT 


The expression for k___ is then 

lih<l 


'‘lei ■ <3 T-'’ 


(el-4388/AT 

(cm ^ in cgs) 

(2.5c^) 

*v 3 — ^ 

= 1.37E-27 G A T ^ n 


, 0.014388/At ,, 

(e -1) 

(m ^ in mks) 

(2.5b) 
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Electron-H Inverse Bremsstrahlung (e + H + hv e + H) 


For the electron-^hydrogen atom absorption coefficient. Ref. 1 

suggests using the results of Stallcop, which are given in terms of the 

electron and hydrogen number densities n„ and n as 

E H 



o o . ^ 

where 0 - h^ m /k, with h^ the ionization energy per unit mass of hydrogen 

I I H 1 

and m^^ the mass of the hydrogen atom. If we put into find 


^LEH 


5/2 

1.43E-45 0^' 
5/2 

1.43E-51 0^ 


2 _T 

T A e ^ n_ n (cm in cgs) (2.7a) 

E H 

1/2 2 _T 

T A e n n (m in mks) (2.7b) 

E H 
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Electron-Neutral Inverse Bremsstrahlung (e + N + hv e + N) 


For other electron-neutral interactions, we use the results of 
John in Ref, 6. .He has provided an expression for the absorption 
coefficient for a large number of neutrals in the form A(T) A per unit 
electron pressure, per neutral particle, including stimulated emission. 

He tabulates a(T) for each neutral. To convert his expression to our k 

Xj 

we multiply by the electron pressure n^ kT and the number density n^ of 
the neutral particle of interest. We also restore the more familiar form 
of the stimulated emission factor by noting that for h^c/AkT small 

j^l - exp (-hpC/XkT)J = h^cA2/kT 


which explains the A^ factor in John’s formula. Then we find, since John’s 
A IS in A, 


k T^A(T) A^(a) (1 - e n 


E N 


LEN 


hp cAA 


Use of cgs or mks units for A and the 


= 9.60 E-1 T^A(T) (1 

LEN 


= 9.60 E-5 T^A(T) A^ (1 


physical constants enables us to write 

-1. 4388/At, , -1 ■ X /o o X 

-e cgs) (2.8a) 

E N 

-0. 014388/At, / -1 . , X /o 

-e )n n^, (m in mks) (2.8b) 

E N 
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The factor A(T) is taken from Table I of Ref. 6. Since it is 
a function of T, we must interpolate. In fact, TA(T) is a smoother 
function of T, so we choose to interpolate in that function. John tabulates 
many neutrals, including Cs, Na, Li, the noble gases, 0 ^, N2, H, 0, N, 
CO, Table 2-1 we give TA(T) for species of present interest. 

Results for 10.6Mm Radiation 


Fig. 2-1 gives a plot of the electron-ion Bremsstrahlung 10.6pm 
absorption coefficient per electron-ion pair, k /n n in cm^, from 
Eqs. (2-5a) and (2-4), vs T. This shows a decrease with T, per electron-ion 
pair. However, the numiber of pairs increases with T, and this dominates 
k as will be seen later. 

XjEiX 

Fig. 2-2 shows the electron-neutral 10.6Pm Bremsstrahlung for H 

5 

and per electron-neutral pair, » from Eqs. (2-7a, 2-6) 

and Eqs. (4-8a) with Table 2-1. These increase with T, but again, their 
behavior will change when the dependence of n^ and n^^ is included. 

Fig. 2-3 gives the electron-neutral 10.6Pm Bremsstrahlung for Cs 

per electron - Cs pair, k__„/n_n„ in cm^, from Eqs. (4 -8a) with Table 2-1. 

JjbN E N 

The minimum reflects the minimum in TA(T) of Table 2-1. The increase of 
electrons with T will change that shape. 
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TABLE 2-1 


TA(T) X 10 


34 


(From Ref . 


6 ) 


Species 

lOOOK 


5000K 

10,000K 

25,000K 

«2 

703 • 


2145, 

2900, 

2800. 

N 

120, 

-1 

500. 

904, 

1748. 

0 

103, 


372. 

654, 

1280. 

CO 

605, 


2250, 

2800. 

3750. 


459, 


1665, 

2930. 

3500. 

°2 

261. 


855. 

1140, 

1900. 

CO2 

1750. 


960, 

1600, 

3250. 

Cs 

57, E3 


24, E3 

28, E3 

30. E3 
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Fig. 2-1 Electron- Ion Bremsstrahlung Absorption Coefficient 
at 10.6 yrn. 
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T (10? K) 

Fig. 2-2 Electron-Neutral Bremsstrahlung Coefficient for at 10.6 Mm. 
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Finally, we inay obtain the actual 10.6ym absorption coefficient 
k by combining the k per pair with the equilibrium composition for a 
H^/Cs mixture as calculated by the method described in Section IV of this 
report. The values from Figs. 2-1, 2-2 and 2-3 are multiplied by the 
appropriate number densities from the equilibritmi calculation and the four 
contributions to k are added to get the total absorption coefficient. The 
result will depend on pressure as well as temperature, because the 
composition depends on pressure. 

When this is done, the results shown in Fig. 2-4 are obtained. The 
solid curves are for a mixture H/Cs = 0,99/0.01, that is, one Cs atom for 
each 99 hydrogen atoms. For comparison, the dashed curves are for pure 
hydrogen. At low temperatures, the large increase in k when cesium is 
added is clearly shown, caused by the presence of electrons from the easily 
ionized cesium. The solid curves flatten when the cesium is fully ionized, 

and even dip slightly as the increasing temperature reduces the density. 

\ 

They rise again when hydrogen ionization begins to supply additional electrons. 
The pure hydrogen absorptxon coefficient is far below that of the seeded 
mixture until near 12,000 K at the highest pressures, and near 9,000 K at 
the lowest pressure. If one consideres 10 ^ cm as a reasonable absorption 
coefficient, then the mixture is absorbing well at temperatures as low as 
2,000 K for 100 atm, 3,000 K for 10 atm, but not untxl 12,000 K for 1 atm. 

The advantage of higher pressure is evident. 

It is important to note that even the seeded mixture does not absorb 
well below 2,000 to 3,000 K. Heating the gas to this temperature must be 
accomplished by some other means, such as introducing a low temperature 
molecular absorber, or by the LSC wave mechanism of conduction and forward 
radiative heating. The cesium is a good absorber of lO.Gym radiation in 
the middle temperature range, and extends the absorption capabilities of 
pure hydrogen down to 2,000 to 3,000 K. 
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2.2 Molecular Absorption 


Laser energy absorption by molecules for the wavelengths of lO.SPm 
and 5.3ym has also been examined by a study of the available literature. 

Since this absorption depends on resonance between an internal molecular 
transition and the laser wavelength of interest, one must match the 
absorber to the laser. Another characteristic of molecular absorption 

is a rapid cut-off with temperature as the molecule dissociates . ^ 

Most molecular absorbers work only at low temperatures, their temperature 
limit depending on the dissociation energy of the molecule. 

We will first consider the absorption coefficients for several 
molecular absorbers, normalized by a fixed moihber density of molecules. 

This will display the intrinsic temperature dependence. Later we will 
conbine this with calculations of the actual number densities in gas 
mixtures, to obtain the total absorption coefficient of the mixture. 

It is important to realize that the molecular absorption coefficients 
discussed below are calculated or measured for conditions of vibrational 
and rotational equilibrium and that they represent averages over a finite 
frequency range (the lines are assumed to overlap or the measurements 
have a relatively coarse resolution) . This could lead to significant 
errors because of the very narrow linewidth of the laser line, particularly 
for absorption by H 2 O where line broadening parameters for broadening 
by are not known, and the line spacing can be very large because of the 
large rotational spacing which can occur. The present work will include 
experimental measurements of absorption of selected 10.6ym laser lines by 
water vapor, made in mixtures at elevated temperatures and pressures. 

The results should help resolve these questions. ^ 
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TUjsorption by H^O at 10>6ym 


Water absorbs lO.SPm radiation, and may be a good low temperature 
absorber, since it doesn't dissociate until about 4,000 K, by which time 

< 

cesium can provide enough electrons for Bremsstrahlung to take over. A 
survey of the literature on this absorption coefficient is shown in Fig. 2-5, 
which gives k vs 1000/T, with T in K and k in cm ^ per amagat, where an 
amagat is 2.69E19 particles per cm^. 

The lines on the right (at low temperature) arise from continuum 

absorption, and are derived from Ref. 7. At 10.6ym, this reference 

2 -1 -1 

presents an absorption coefficient for 296K as 2.25 E-22 cm atm molecule , 
and gives a temperature correction factor as exp [l800 (1/T - 1/296 . 

For an amagat of particles, this becomes 

k = 1.38 E-5 P (atm) (cm"^ ama 

H2O 

which shows a pressure dependence, unlike the line absorption mechanisms, 
whose only pressure dependence is in the particle density. If one converts 

to cm ^ by using a water number density of n , then the absorption 

. . 2 
coefficient becomes 


= 5.14 

E-25 

1800/T n 


(cm“^) 

- 3.77 

E-3 

ieoo/T,„ J 

( H.,0) 

/ T 

( 

= 7.00 

E-47 

V 2 / 

1800/T / \2 

" V"h,o) 

T 
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Fig. 2-5 Absorption Coefficient of HO at 10.6 Hm. 



The squeure dependence on pressure or number density for this absorption 
mechanism indicates an advantage in going to higher pressure for low 
temperature absorption of 10.6ym radiation by water. 

The left-hand (high temperature) dashed curve is from Ref. 8, 
and is due to vibration-rotation line absorption. This curve is identical 
to that of Ref. 9, which was probably calculated the same way. Recently 
Fowler, Newman and Smith, in Ref. 10, have also calculated this contri- 
bution for water, and the ic curve in their Pig. 23 is quite close to the 
dashed curve of Fig. 2-5. In fact, the points in their Fig. 23 referenced 
as Ludwig's data are actually calculations from the same source as the 
dashed curve in Fig. 2-5. 

There are also rotational lines whose wings make a contribution 
at 10.6ym. They have been calculated in Ref. 11, but give very low 
values near 600K. 

The few measurements of water vapor absorption in the literature 
are also plotted in Fig. 2-5. There is a group at low temperature 
between 400 and 500K from Refs. 12 and 13, and a group between 1535 and 
2200K. These latter are calculated from spectral emissivities measured 
in Ref. 14. They fall very close to the dashed curve. All these measure- 
ments have coarse resolution in frequency. 

The results of a calculation based on the "McClatchy tape" au:e 
given in Ref. 15, and they have higher valuhs of k at low temperature. 
However, not enough details are given there to inspire confidence in the 
results. AtPSI, a preliminary analysis performed by B. David Green, also 
using the line-by-line information in the McClatchy tape, indicated low 
absorption coefficients near 600K. 
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Based on all this information, we have drawn an upper-bound 
estimate for the line absorption, shown as the chain-dashed curve in Fig. 2-5, 
which was the curve used for our calculations, when added to the continutam 
absorption- 

From this study, we conclude that water vapor absorption of lO.Gym 
radiation has a minimiam around 500K, and is not as high as would be 
desirable below lOOOK. However, the measurements are sparse and scattered, 
have coarse frequency resolution, and none have been made between 1500K and 
500K, where the minimum is indicated in Fig. 2-5. It is clear that some 
careful, high resolution measurements would be very valuable in determining 
whether water vapor is a useful low temperature absorber for 10. Gym 
radiation. 

If water vapor proves to be an unsatisfactory absorber at low 
temperatures, there are a nxirober of other molecules which can be considered. 
Several have been studied theoretically in Ref. 10, and experimentally 
with results not yet reported. But in addition to the absorption coefficient, 
one must know the chemical stability of these molecules before their true 
potential as absorbers can be evaluated. This requires knowledge of their 
kinetics as well as their equilibrixam properties with hydrogen. Study of 
these other 10. Gym abosrbers is beyond the scope of the present work. 
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Absorption by H 2 O at 5.3Pm 


The absorption coefficient of H.O for 5.31Jin radiation is shown 

-1 -1 ^ 

in Fig. 2-6, where k in cm amagat is plotted against 1000/T. This 
curve is taken from Ref. 16, up to 3000K (1000/T = 1/3). For higher 
values of T, the curve is extrapolated horizontally. It is seen that water 
is a good absorber over the whole range from 300 to 3000K, since k varies 
only from 0.06 to 0.22. 

Absorption by NO at 5.3ym 

The absorption coefficient of NO for 5.3Vm radiation is shown in 
Fig. 2-7, where k in cm ^ amagat ^ is plotted against 1000/T. This 
curve is taken from Ref. 16 up to 3000K, and extrapolated horizontally 
for higher values of T. The shape of this curve shows excellent absorption, 
above 1.0, at low temperatures, and near 0.5 at high temperatures. NO is 
not a thermodynamically stable species lander almost any conditions of 
interest; however, below about 2500K the chemical kinetic rates for its 
destruction are so slow as to make it a candidate for a low temperature 
absorber. In the scenario of interest it will form H 2 O and N 2 when it 
dissociates so that there would be a good absorption coefficient from 
300K to 3000K if NO is initially present. 

Absorption by CO at 5.31-im 

The absorption coefficient of CO for 5.3)Jm radiation is shown in 
Fig. 2-8, where k in cm ^ amagat ^ is plotted against 1000/T. This curve 
is taken from Refs. 16 and 17. It covers the temperature range from 
1000 to 6000K, over which k varies from 2.3E-3 to 0.31. The agreement 
between the values provided by the two different references gives confidence 
in the prediction of k for CO. Fig. 2-8 shows that only above 2000K may 
CO be considered a good absorber of 5.3Mm radiation, so it acts as a high 
temperature absorber, while NO and H 2 O are low temperature absorbers. CO 
is stable to very high temperatures (> 6000K) . 
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Absorption Coefficient of NO at 5,3Pm 
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Results for Gas Mixtures 


The absorption coefficients for individual species can be combined 
with chemical calculations of the composition of a gas mixture, ,to provide 
the total absorption coefficient for the mixture at a given pressure and 
teirperature . We need only multiply each absorption coefficient by the 
nvnnber density of the one or two species which partake in the process 
leading to that absorption mechanism. Then all the coefficients are 
added to give the mixture absorption coefficient k . 

The composition of the gas mixture was calculated by an existing 
PSI computer program which accepts a list of species, a pressure and a 
temperatTire and calculates the equilibrium composition, enthalpy and density. 
This program can also perform calculations with some species frozen, by 
omitting from ,the species list those species which would result from 
reaction of the species it is desired to freeze. 

Calculations have been made for both 10.6Pm and 5.3Vim radiation 
at pressures of 1, 3, 10, 30 and 100 atm, for temperatures up to 6000K. 

For 10.6Pm, the mixture considered was H 2 /H 2 O/CS in the six composition 
given by 

H^/H^O/Cs = 0.945/0.05/0.005, 0.940/0.05/0.01 

= 0.965/0.03/0.005, 0.960/0.03/0.01 

= 0.985/0.01/0.005, 0.980/0.01/0.01 

The absorption coefficients are shown for 3, 30 and 100 atm in 
Figs. 2-9 and 2-10 for all six compositions. At 3 atm, the absorption 
coefficient never rises above 0.1 cm ^ with these compositions and is below 
10 ^ cm ^ for temperatures lower than lOOOK. The solid curves have 0.005 Cs 
while the dased cturves have 0.01 Cs. The extra cesium makes a noticeable 
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increase in k above 2000K, where the cesitmi begins to produce electrons 
L 

SO Bremsstrahlung can be effective. However, below that temperature, all 
the absorption comes from water, and both sets of curves are the same 
because they have the same amount of water. 

The effect of pressure in increasing k is shown by Fig, 2-10, 

Xj 

where the pressure is 30 and 100 atm, rather than 3 atm. Here the 

highest values of k are near 25 cm ^ for the 0.01 Cs mixtures and 14 cm ^ 

for the 0.005 Cs mixture. These represent factors of 500 to 700 increase 

over the corresponding values at 3 atm and show the very favorable increases 

in k that are obtained with the factor of 33 increase in pressure. 

L 

Again, the pure water absorption is the sole contributor to k below 

Li 

2000K. 


In Fig. 2-10, the minimum of water absorption around 500K is clearly 
seen, with values of 0.002 to 0.02 at 100 atm for these mixtures. This 

present some problems for absorption, so it is important to verify that 
the minimum shown on Fig, 2.5 actually is correct* If it proves correct, 
then other absorbers for lO.GMm should be studied, since absorption 
coefficients below 0.01 are probably not desirable. 

At the high temperature end, the 100 atm absorption coefficients 
are more than adequate; it is probably not necessary to operate at such a 
high pressure. The 30 atm curves show maximum values of k of 1 to 3 cm 

Li 

From these results we can conclude that a 0.005 or 0.01 Cs mixture 
with a few percent H 2 O in H^ provides good absorption of lO.GPm radiation 
at 30 atm above 2500K. However, below that temperature, the absorption 
falls rapidly if our present estimates of water absorption are correct, and 
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better absorption is needed than that provided by up to 0.05 H 2 O. Either 
other absorbers should be studied, or considerably more water is needed. 
The addition of water (and consequent subtraction of H 2 ) increases the 
molecular weight of the mixture and reduces the specific impulse which 
can be attained. Therefore, there is a tradeoff between absorption 
coefficient and amount of molecular absorber. 

The calculations for 5.3]Jm radiation considered two mixtures. 

Both used CO as the high temperature absorber, but one used water as the 
low temperature absorber and the other used NO. The mixtures considered 
were 


H2/CO/H2O = 


0.90/0.05/0.05 

0,94/0.03/0.03 

0.98/0.01/0.01 


H /CO/NO 
2 


= 0.90/0.05/0.05 

= 0.94/0.03/0.03 

= 0.98/0.01/0.01 


For each mixture, two calculations were made. In the first, all species were 
assumed to be in equilibrium. In the second, the NO or the H 2 O was taken as 
frozen, unable to react. These provide the two limits between which the 
actual chemical state will vary. It is important to consider these limits 
because NO is not a true equilibrium state at low temperatures, and neither 
is H 2 O in the presence of CO. However, both NO and N 2 O will actually react 
very slowly, and in laser heating applications, we do not expect true equi- 
librium to prevail. 
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The equilibrixom composition calculations show this lack of H^O 
at temperature below lOOOK. The water is replaced by CO^ in this region. 

When we put in NO instead of water, the composition is nearly the same. 

The NO is replaced by N^ over the whole tenperature range, and the O 
becomes water above lOOOK, and CO^ below there. As far as the equilibrixim 
species go, the H^/CO/NO and N^O/CO/H^O mixtures are very similar. 

In both of them, the CO is present at a constant mole fraction above lOOOK, 
and so can serve as the high temperature absorber. 

When the composition is combined with the individual species 
absorption coefficients we get the k curves shown in Figs. 2-11 (3 atm) 

Jj 

and 2-12 (100 atm) . 

The results for 3 atm show values of k mostly between 0,01 and 

-4 -1 ^ 

10 cm . The pure equilibrium cases are nearly the same whether NO or 

H^O is used. They fall off sharply below lOOOK because of the disappearance 

of H^O. The frozen water curves hold up better below lOOOK, though they 

also fall off below 500K because of the dip in water absorption. They are 

also slightly higher at the high temperature end. The frozen NO curves 

are somewhat higher everywhere than the frozen water curves, and they rise 

sharply at low temperature instead of decaying. This reflects the high 

absorption of NO at low temperatures . 

The results for 100 atm in Fig. 2-12 are qualitatively very similar 

to the 3 atm case, but all values are about a factor of 30 higher, 

corresponding to the pressure increase. This improvement is not as great 

as for the lO.OPm case, because the molecular absorption is proportional 

only to the first power of number density, while for 10.6pm, the Bremsstrahlung 

is proportional to the product of the electron and molecular number densities . 

The values of k for 100 atm at the higher seed fractions are mostly between 
- 1 ^ 

1 and 0.01 cm , which are adequate absorption coefficients, though not as 
high as for lO.Gym at the same pressure. 
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The high values of k at low temperature provided by frozen NO 

L 

would be very useful if the chemical kinetics will permit this composition. 

For 5.3Vim radiation we conclude that frozen NO will permit good 
absorption at low temperature; frozen water will give lower, though 
probably adequate absorption. The pressure level must be higher to reach 
the same values of k as for 10 ,6ym absorption. The trend of absorption 

Xj 

with increasing temperature is down for 5.3ym, since no new absorbing 
species are being created, while the density is decreasing. This is in 
contrast to lO.GPm, where the trend is up because of the creation of 
electrons . 
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III. RADIATIVE AND ABSORPTIVE PROPERTIES MODELING 




The use of CW laser radiation to heat gases for propulsion is 
advantageous because high temperatures can be reached, leading to the possibility 
of high specific impulse. However, the high temperatures also lead to sig- 
nificant radiation from the gases, which must be handled by the cooling system 
of the thruster. In the continuous heating mode, where molecular absorbers 
are used, temperatures up to 6000K can be reached, at pressure levels of 
up to 100 atm. At temperatures above 2000K, radiation from the vibrational 
modes of the molecular species in the gas will be the important contributor 
to the total radiation. For hydrogen, or an alkali metal seed, such as cesium, 
which ionizes, the free-free (Bremsstrahlung) and free-bound radiation will 
also be significant at higher temperatures. 

In the LSC wave heating mode, where much higher temperatures 
are reached, line radiation may also be important. 

The molecular species considered in this study include NO, CO 
and H^O, as well as the primary propellant At the higher temperatures 

OH is also formed and must be considered. When NO is used as an absorber, 
finite rate kinetic calculations show that it will react with H 2 to form 
H 2 O and N 2 at temperatures above 2000K, so radiation from NO is not important 
and has not been calculated. However, such a calculation could also be made 
by the procedure described below for the other radiators considered. The 
species CO 2 is not formed in significant quantities in the mixtures studied 
so far, but its radiation was also calculated for two reasons. First, it 
might be a useful absorber for future consideration, and second, there are 
data and calculations available for CO 2 radiation which can be compared with 
the results of the calculation method used here. 
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The radiation from the species U^, H, H , E, Cs, Cs^, CO, 

^^2' ^2^ have been calculated here, including molecular radiation, 

free- free and free-bound contributions. The method used for the molecular 
radiation is that described in the books by Penner (Ref. 18) and Penner 
and Olfe (Ref. 19). It is called the "just-overlapping line" or "continuously- 
overlapping-line" approximation. This method is particularly applicable 
at high pressure where lines are broad, for large optical path where the 
wings of lines become more important, or at high temperature where the number 
of significant lines in a wavelength interval increases. The purpose of 
this work is to develop a radiation model which is simple enough to be easily 
implemented on a computer to calculate the emission from a mixture of the 
gas species considered, both spectral and total. The results can then be 
used to construct even simpler and less detailed models of emission for use 
in the thruster and LSC wave computer programs. We do not require very high 
accuracy for the models, but desire to obtain the emission to what might 
be called engineering accuracy. 

The radiation from many of these species has been studied before, 
of course. For the combustion gases CO, CO^ and H^O, Hottel long ago produced 
his famous radiation charts (Ref. 20) based on experiments. However, they 
were directed to furnace applications, and were limited to temperatures below 
2000K and pressure-length products below 10 atm- feet. We are interested 
in temperatures up to 6000K and higher values of the pressure-length product. 

For these species, then, our work may be looked on as an extension of the 
Hottel charts to higher temperatures and pressures. Pure hydrogen is also 
a very well-studied gas, and up to 10,000K the paper of Olfe (Ref. 21) 
is particularly useful. 

We will see that the results of our model compare very well with 
the calculations and data of other investigators, and give confidence in 
the present method of calculating molecular vibrational radiation up to at 
least 6000K. 
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3.1 Molecular Vibrational Radiation 


Band Strenths for Vibrators 


The band strength of a vibrational band (often called the integrated 

absorption) is an important distinguishing characteristic of the band. We shall 

-1 -2 

use the units amagat cm where an amagat is the number density at 1 atm 

"3 

and 273. 16K, which is 2.687E19 cm . Band strengths are normally given at 
room temperature and their temperature dependence must be determined. It 
is readily found by use of a simple harmonic oscillator approximation for 
the vibrator, which will now be described. 


When a vibrator makes a transition from a lower state Z specified 
by the quantum number n^, n^f-** to an upper state u specified by n^ + 6^, 

^2 band strength is 


L c 

a (n. n. + 6.) = — 

^ ^ 8ttv2 


-E„/kT 






(3-1) 


Here c is the speed of light, is Loschmidt's nvunber, is the energy level 
of the lower state relative to the ground state, hv = energy 

change in the transition g is the degeneracy of the upper state, Q is the 

V 

vibrational partition function of the molecule and is the Einstein 
coefficient for spontaneous emission for the transition. 

in the harmonic oscillator approximation, the populations, degeneracies 
and Einstein coefficients of the higher level are related to those of the 
lower level by simple rules so that Eq. (3-1) can be summed over the harmonic 
oscillator states to give the total band intensity in terms of the properties 
of the lowest absorbing level. The result for a band is 
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a = \ (n . -> n. + 6.) 
T /_j 111 

n . 

1 


8ttv 


A g f(T) 
uJl u 


f (T) 


r, /, m A-i -hpV/kT, 
^-E^/kT (1-e ^ ) 


n (l-e 


(3-2) 


Here v = where cu)^ is the vibrational frequency of mode i, so that 

(jO^ is the corresponding wave number. 

In terms of the total band intensity at 300K, that at any temperature 
can be written as 

cx^ (T) = (300) f(T)/f(300) (3-3) 


At room temperature hpV/kT is usually large compared to unity, and for transitions 
from the ground state of the molecule f(300)^ 1. 


) 
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Emissivity 


As radiation traverses a gas, some is re-absorbed, so the radiation 
emitted from the end of a column of gas is less than that for emission only* 

The absorption depends on the optical path and the absorption coefficient. 

We here define the optical path length x the physical length times the 
number density of the radiator in amagats, and define the absorption coefficient 
P so that PX is dimensionless. 


The radiation emitted from the end of a homogeneous column of 

-PX 

gas is the blackbody function times the absorption factor 1 -e . This is 

4 

usually expressed as an emissivity e times the blackbody flux OT , where 0 
is the Stefan-Boltzmann constant. The emissivity is used as the engineering 
expression for the emitted radiation. 


£ 




(1-e 



(3-4) 


Here the blackbody function is defined in wave number (U3 = v/c) space as 


B 


to 


27ThpC ^0)^ 


(3-5) 
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4 

such that its integral over U) from 0 to ~ is OT . In order to find e we must 
evaluate the integral in Eq. (3-4) with the proper variation of P with respect 
to CO. 


We will evaluate the integral using the continuous-overlapping- 
line approximation. First we apply Eq. (3.4) to a band, centering the integral 
on the wave number co of the band center. We take the blackbody function out- 
side the integral, evaluating it at the band center because the band width 
is small compared to the variation of B^. The band is approximated as symmetric 
about CO, and the integral is extended to “ since the contributions in the 
band wings dies off rapidly. The result is 


e 



I 


(l-e d(o)-a)) 


(3-6) 


The absorption P is taken to be the average of the line strengths 

for corresponding lines in the R and P branches i divided by the line spacing 

2b , and so becomes 
e 


P 


2j +1 
4 Be 




-Yj(D+l) 


(3-7a) 


Here a„ 
T 


is the band intensity, y is defined by 


'Y = h^cB /kT 
P e 


(3-7b) 
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and j is the rotational quantvim nvunber of the upper state. 


We treat j as a continuous variable, related to o) by the line 
spacing U) = O) + 2jB^ and introduce the variable 

^ = Y j (j+1) + Y/4 = Y(j+i)2^ 


so that dC =Vy dj, 


dj = d {(A)- 0))/2B . 

e 


Then 


P 


Y/4 - 5 " 

2B e e 


In terms of C the expression (3-6 ) for e can be written 



(3-8) 


The one parameter integral I (K) can be evaluated numerically. 
As K 0, I ■+■ K. For K large, an asymptotic analysis shows that I goes 
like (S-n K)^. A good fit to the integral is 
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I 


7.2 [1- exp ((-K/7.2)], K < 20 


= 4.13 


[iln 0.721K] 


i 


K > 20 


(3-9) 


With this approximation for I (K) , 

in terms of T, B , a and X by Eqs. (3-8) , for 
e T 

are found from Eqs. (3-2) and (3-3), while 
to the physical path length I by 


the emissivity is expressed 
each band. The band strengths 
the optical path X is related 


X = £. [273.16 p(atm)/T] 


(3-10) 


The factor in square brackets is the nimber of amagats at the desired pressure 
and temperature. Once the bands and other physical constants of the molecule 
of interest are described, the emissivity can be easily calculated. It is 
a function of T and the pressure-length product pS,. 

For molecules with multiple bands, the emissivities are added 
unless the bands overlap. In the case of overlapping bands, a reasonable 
approximation is to use the sum of the overlapping band strengths. 

In flow problems such as the laser heated thruster, radiation 
models are needed to describe radiant emission, and possibly also radiation 
transport. To construct such models, one needs to know the spectral emissivity 
e^, which is the integrand of Eq. (3-4) . That can be calculated from the 
relations already given. The blackbody function is given in Eq. (3-5) , the 
band strength in Eqs. (3-2) and (3-3), and the absorption P in Eqs. (3-7). 

The rotational quantum number j is related to the wave number and the band 
center wave number w by 
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w — w 


(3-11) 


= 2jB 

e 

It should be recalled that X is length times the number density of the radiator 
in amagats, so that P has that number density built in. If one wants the 
absorption coefficient for a number density n of radiators, it is found from 


k = nP/2.6868E19 


where k is in cm ^ and n is in cm 


(3-12) 


Emissivity of CO and OH 

For diatomic molecules there is only one vibrational mode, and 
a small number of bands. For OH we consider only the fundamental band, while 
for CO we consider the fundamental and the first and second overtone bands. 
For both molecules, the lower state energy E^ = 0. The parameters for OH 
and CO are given in Table 3-1. 


The calculated emissivities of OH and CO are given in Figs. 3-1 
and 3-2 as functions of T for various values of pi.. The dependence on pA 
is nearly linear for OH but shows much less than linear dependence for CO, 
reflecting the strong effect of reabsorption for this molecule. 
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TABLE 3-1 


Parameters for OH and CO 


Molecule 

o)(cm 

, -lx 
B (cm ) 

6 

-1 -2 
a (300) (ama cm ) 



e 



OH 

3570 

18 

1 

100 

CO 

2143 

1.9 

1 

260 


2 2.2 

3 0.011 

TABLE 3-2 

Parameters for CO^ 

cm ^ ^3 ” 2350 cm 

a(3Q0) (ama ^ cm 

2770 
240 
67 
7.5 
1.9 
0.16 

50 


B =0 .39 cm 



0 

0 

1 

1 

0 

0 


0)^ = 1388 cm 


' 0 1 

1 0 

0 1 

-1 0 

4 1 

3 0 




Fig. 3-1. Calculated Emissivity of OH 
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Emissivity of CO 2 


For the triatomic molecule CO^i the calculation of emissivity 
is significantly more complicated than for the diatomics. The parameters 
used are given in Table 3-2. The frequency of each band is given by 

V = c '5 . w. 

^11 

The band with the negative value of 6^ corresponds to a series of transitions, 
the lowest of which is from the 100 state to the 010 state, and has a strong 
temperature dependence of the band intensity a. 

In Eq. (2.4), E^^ corresponds to the state given by the negative 

values of 6; 

E . = - he 6 . w . 

ii ^ Zu ^ ^ 

6<0 

If there are no negative values of 6 then E^^ = 0 . 

Another property of CO 2 which must be taken into account is the 
natural abundance of the carbon 13 isotope, which rs about 1%. For CO, the 
data on emissivity is not certain enough and the bands not strong enough to 
make this worth considering, but for CO 2 the effect is considerable. The 
bands of isotopic CO are shifted enough from normal CO_ that they do not 
overlap for conditions of interest here, but they are in the same spectral 
region. The effect of the isotope is therefore approximated by assuming the 
same band intensities, but at 1% of the partial pressure. We therefore add 
to the emissivity of normal CO 2 at pil the emissivity of CO 2 at p£/100. 


53 




At large pZ where the CO^ radiation becomes highly absorbed, this becomes 
a very important contribution to the emissivity, though always less than a 
factor of two increase. 

The calculated emissivities of CO^ are shown in Fig. 3-3. They 
decrease with temperature, except for a slight maximum at the lowest value 
of p£- 


Emissivity of H^O 


The calculation of the emissivity of H^O is complicated by the 
difficulty in describing the rotational structure of the bands. We have taken 
the width to be given by the variation of intensity with the principal quantum 
number. The band envelope then approximates that of a linear molecule but 
with the rotational constant given by the square root of the product of the 
two layer rotational constants of H^O. We have also combined the two vibrational 
stretch modes of H^O into one to simplify the calculation, leaving only two 
modes. The pure rotational band was omitted because it contributes significantly 
to the emissivity only below the temperature range where radiation is import- 
ant for our application. 

The values of the parameters for H^O used in the calculations 
are given in Table 3-3. There are five bands, each with two transitions. 

The calculated emissivities are given in Fig. 3-4. The off 

below lOOOK is partly due to our neglect of the pure rotational contributions. 
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TABLE 3-3 


Parameters for H 2 O 


B = 11.8 cm 

e 


0 

1 

1 

2 

2 


= 3655 cm 



1 

0 

1 

0 

1 


= 1596 cm 

a(300) (ama ^ cm 

300 

200 

24 

18 

1.0 


TABLE 3-4 

Vibrational Emissivjty of Molecular Hydrogen at 
a Total Pressure of 100 atm and Path Length of 30 cm 


T(k) 


n (cm 

“2 


p (atm) 

"2 


£ (Eg. 3-14) 


e (Ref. 3-14) 


1000 

7.34 E20 

100. 

1.14 E-2 

1.00 E-2 

2000 

3.67 E20 

100. 

5.83 E-3 

5.34 E-3 

3000 

2.41 E20 

98.4 

1.90 E-3 

1.92 E-3 

4000 

1.56 E20 

85.3 

5.55 E-4 

5.59 E-4 

5000 

7.78 E19 

53.0 

9.66 E-5 

1.56 E-4 
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Emissivity of 


The vibrational radiation of the H 2 molecule has been considered 
by Olfe in Ref. 21. He has calculated the band intensity. From his Fig. 5 
and Eqs. (10) we can write 

a = 0.019T~°‘^ P„ (atm) (3-13) 

T H2 

The parameters for are w = 4200 cm”^, = 60 cm~^. 

If we use these values in Eqs. (3-8) , and (3-10) we can 
calculate e for H 2 vibration from the expressions 

_ 3.16E11 I (K) 

6043 /T , 3.5 (3-14) 

e -1 T 

K = 0.8037 Ip^ /T^‘^ 

H2 

Olfe presents a calculation for a total hydrogen pressure of 100 atm, 
5. = 30 cm, and identifies the vibrational contribution in his Fig. 16. By 
calculating the equilibrium state of hydrogen at this pressure and any T we 
can find the number density of H 2 , and so its partial pressure. Then 
Eq. (3-14) permits a calculation of e for ij, = 30 cm. The results are given 
in Table 3-4, together with the values for vibration from Olfe's Fig. 16. 

The agreement is quite satisfactory up to 4000K, and above that temperature 
the vibrational contribution is unimportant. So Eq. (3.14) can be used for 
the total emissivity of hydrogen vibration. 
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Eqs. 

is 


The absorption coefficient is found from Eq. (3-12) , using 

(3-7a) , (3-11) , and (3-13) and the relation between p and n , which 

”2 H2 


p = 1.3625E-22 n„ T 


(3-15) 


when p is in atm. 

“2 


Comparison with Other Calculations and Experiments 

There are other calculations and also experimental data with which 
to compare the results of the present model of molecular emissivities. 

For CO, Fig. 3-5 shows the limited data of Hottel (Ref. 20) and 
the radiation in the thin limit obtained from the calculations of Young (Ref. 
22), as well as the present calculations from Fig. 3-2. Young has performed 
detailed line-by-line calculations. Our results for small optical path are 
with 10% of Young's over the entire temperature range considered. This in 
effect justifies, at least for CO, the simple harmonic oscillator approximation 
for obtaining total emissivity. 

Hottel 's data is in good agreement with our calculations at 
pi. = 1 atm-ft. and fairly good agreement at 0.1 atm-ft. For 0.01 at ft., 
the data must be regarded with considerable suspicion, since it is in 
sharp disagreement with both the present results and those of Young. 

Based on the comparisions of Fig. 3—5, the present model for CO 
emissivity is considered satisfactory. 
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Fig. 3 - 5 . Emissivity of GO 
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Figure 3-6 presents the data of Hottel (Ref. 20) and the cal- 
culations of Fenner and Varanasi (Ref. 23) together with the present calculations 
for C02* The agreement is very good, and again, gives confidence in the present 
model. Hottel gave correction factors for various pressures which account 
for non-overlapping line effects. Our calculations therefore correspond to 
the high pressure, continuously-overlapping line results. 

Data and calculations for H2O are presented in Fig. 3-7. The 
calculations of Fenner and Varanasi (Ref. 24) include the effects of finite 
line width in a calculation very similar to ours. The Ludwig and Ferriso 
upper limit points are quoted in Ref. 25 as coming from Ref. 26. The Ludwig 
and Ferriso points are from their calculations of the integral over frequency 
of results from a band model in Ref. 25. Both these latter results also 
appear in Ref. 27. Figure 3-7 shows our present calculations to be in good 
agreement with the other work, particularly at the higher temperatures and/or 
lower pi. The main discrepancy can be attributed to the difficulty in character- 
izing the band width. This problem decreases with increasing temperature 
because lines associated with higher values of the secondary rotational quantum 
number fill in between the principal lines at higher temperature, thus making 
the overlapping line approximation much better. 
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3.2 Electronic Radiation Processes 


The reverse of the processes which produce absorption in the presence 
of electrons will produce radiation. These include two free-free processes: 
one between neutrals and electrons (electron-neutral Bremsstrahlung) and one 
between ions and electrons (electron-ion Bremsstrahlung) . In addition, there 
are the bound-free processes which include photo-detachment and photo- 
ionization. We will consider these for both hydrogen, which is the major 
gas species, and cesium which is the alkali metal seed being used in the present 
study. 


The radiation will be formulated in terms of absorption coefficients 
for each process, which is a function of the temperature, composition, 
wave number, and physical properties of the species. From them, the spectral 
emissivity and total emissivity e can be obtained for a given path length 
by means of the relations 


e 


= / e doj, z 

J 03 0) 


B 

03 



OT 


4 


(3-16) 


where B^ is the blackbody function already defined in Eq. (3-5) . This is the 
emissivity normal to the center of the base of a hemisphere of radius 

The absorption coefficients, in turn, are of the form of a cross- 
section times the number densities of the species participating in the process 
times the stimulated emission factor 1 - exp (-hpC03/kT) . Below we will present 
the cross-sections (denoted by o) for the various processes of importance. 


/ 
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“PAGE MISSING FROM AVAILABLE VERSION” 

C5'U 



For hydrogen molecules, Table 2-1 shows that T A(T) may be 

1 45 

from 1000 to 25000K by 3.7 E-33T ' so that the cross-section is 


agH2 = 3.54E-33T^*'^^ 


and the absorption coefficient is 


-3 1 45 

k ^ « 3.54E-33 (a3 
EH2 ^ 4. 


( - -h_ c(A)/kT I / -l\ 

l-e P ) 


fitted 


(3- 19b) 
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Free-Bound Radiation 


Free-bound radiation is caused by free electrons being captured 
by an ion or neutral atom to form a neutral atom or negative ion with the 
release of energy as radiation. The absorption coefficient is readily found 
for the bound- free process for neutral hydrogen and cesium, and then used 
to find the emissivity for the free-bound (reverse) process. 


For hydrogen atoms, the bound-free cross-section is given by the 
sum of the cross-sections for the available bound states 


a = 7.9E-18 

H 


Z (s) --Aw- 

>/ryM 


N>/Ry/a3' 

~N 

max 

z 

N = 1 



where Ry is the Rydberg (109737 cm and is the principal quantum number 

of the electronic state of hydrogen which will overlap that of an adjacent 

atom. For simplicity we took N = 40. 

max 

The absorption coefficient for free-bound radiation of hydrogen 

atoms is found from 0 by using their number density n in the expression 

H H 


, ^ /■■ «»“h^ca)/kT. , -1, 

(1 -e p ) (cm ) 


(3-20) 


For bound-free absorption of cesium, the cross-section can be ap- 
proximated adequately by the quantum defect method, which we treat as a perturbation. 
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The quantum defect method defines the energy level of an electronic state 

below the ionization level in terms of a modified principal quantum number. 

* 

This modified quantum number is denoted by n and is found to be related 
to the principal quantum number by a defect <5^^^ as 

n* . n - 

Where is only a function of the angular momentum quantum number, Z. For 
cesium values of 6^ are given in Table 3-5, taken from Ref. 29. 

The energy for ionization from an electronic state is then given 
by 1 


V ■ 

The ionization from the ground electronic state (n = 6, 1!, = 0) is then 

Wgo = Ry/(1.87)^ = 0.287 Ry , 

and the electronic partition function becomes 

5 "^maXjA r he *1 

Qe£ “23 ■ 23 2(2JI + 1) exp - — (0.286 Ry - 

Z=o > n^ L J 
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where n„ is also given in Table 3.5 as is the upper limit on n, n .. 

itidx g Xf 

The bound-free cross-section for each electronic state, (n,Jl) 
is then given in the same way as for hydrogen (we performed some reduction 
of the expression) , and is given by 


a 


n,Z 


7.9 E-18 



E 7.9 E-18 





except that for the ground electronic state (n-=6, 5,= 0) we do not allow 
a (6,0) to be greater than 2 E-18. This then allows acceptable agreement 

with the limited data of Allen, Ref. 30. 


The total cross-section for bound-free absorption by neutral cesium 


is then 


Cs 


5 max, £ 

i: 2 

£=o n> n. 


2(2Jl + 1) a . exp [- r— (0,286 Ry - to 

n , A/ K.T 


to > (0 


n, £ 


and the absorption coefficient is related to the cesium number density n^^ 
by 


Cs 


/-i -h cto/kT. , -1. 

= o„ n„ (1 -e P ) (cm ) 

Cs Cs 


(3-21) 
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The bound- free absorption by negative hydrogen ions is importcint 
when cesium is present, because electrons provided by the cesium attach 
to the hydrogen atoms. The photo-detachment cross-section for this pro- 
cess is given in Ref. 31. We used this data as a table, given in T2QDle 

3-6. The cross-section a was interpolated linearly in this table, 

HM 

using log^ _ 0) and a as the variables. Then the absorption coefficient 
10 HM 

is related to the number density n^ of negative hydrogen ions by 

-h co)/kT \ 

) (cm"^) (3-22) 


^HM ^HM 


”hm 
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TABLE 3-6 


Cross-Section for Photo-Detachment of H (Ref. 31) 
U) (cm X 10^^ (cm^) 

HM 


1000 

0 

6050 

0 

6250 

0.3 

6670 

1.0 

7140 

1.59 

8330 

3.4 

9090 

4.0 

10000 

4.3 

11110 

4.4 

12500 

4.4 

14300 

4.3 

16670 

4.0 

20000 

3.3 

25000 

2.7 

33330 

1.9 

50000 

1.1 

10000 

0.3 
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Results 


1 


1 


1 'ij w »-i .} ! 


With the absorption coefficients known for the various electronic 
radiation processes the spectral and total emissivities can be calculated 
from Eq. (3-16). 

For pure hydrogen, we have considered six different processes; 
electron- ion, electron-H, and electron-H^Bremsstrahlung, H and H free- 
bound and molecular vibration, j The correspondjing absorption coefficients 
are given by Eqs. (3-17, 3-18, 3-19b, 3-20 and 3-22) for Bremsstrahlung 
and by Eqs. (3-12, 3-7a, 3-11, 3-13 and 3-15) for vibration. The six 
absorption coefficients are added to obtain Zk^ for use in Eq. (3-15) . 

The total emissivity e for 100 atm, and £ = 30 cm is shown in 

Fig. 3-8 for T up to lOOOOK as the dashed curve. Also shown as the solid 

curve are the results of Olfe (Ref. 21, Fig. 16). The agreement is very 

good down to 2000K. Below there, Olfe's results are higher because he included 

a rotational contribution which we have neglected. However, the radiation 

4 2 

below 2000K is very small because a T =5.67 W/cm at lOOOK, so our error 
there is unimportant. Therefore, the comparison shown in Pig. 3-8 indicates 
our calculation for hydrogen radiation is quite satisfactory up to at least 
lOOOOK. ■ 


For pure cesium, we have three contributions to the absorption 
coefficient: electron-ion and electron-Cs Bremsstrahlung and Cs free-bound. 
They are given by Eqs. (3-17, 3-19a, 3-21). These three are added to give 
Zk^ in Eq. (3-15). To assess the importance of cesium as a radiator, we 
first plot the normalized blackbody function T in Fig. 3-9 as a function 

of 0), for various values of T. Then, in Fig. 3-10, we show the absorption 
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.g. 3-8. Emissivity of Hydrogen forp=100 atm, 5,=30 cm. 
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Fig • 3“9 . Normalized Blackbody Spectral Function 
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Transparent spectral Emissivity of Cesium for 
p = 1 atm, fi, = 1 cm. 
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4 

coefficient k = Zk. times B /a T on the same scale for T = 5000, 7000 and 

1 0 ) 

lOOOOK. This is the spectral emissivity for Z small enough so that kil is 
small compared to unity. These calculations are for 1 atm, which is the 
highest partial pressure of Cs we expect to deal with, since it corresponds 
to 1% cesixam in 100 atm of hydrogen. 

Comparison of Figs. 3-9 and 3-10 can give an indication of whether 
the radiation for any path length Z at any given wave number is black or trans- 
parent, or between those two limits. If the value on Fig. 3-10 exceeds the 
corresponding value on Fig. 3-9, then for = 1 cm, k£ > 1, and the 

radiation is close to black. For example, this is the case for 5000K at 

2E4 cm On the other hand, at lOOOOK, 1E4 cm ^ the value on Fig. 3-10 
is much smaller than that on Fig. 3-9, showing that k£<<l for Z = 1 cm. 

If £ = 10 cm, all the values on Fig. 3-10 must be raised by a factor 

4 

of 10 to yield k£B^/aT . Then most of the spectrum has k£>l, and the radiation 
will be mostly blackbody. 

For the approximately cylindrical geometry of the thruster; the 
path length of importance is the radius (if the cylinder is not short and 

fat) . This is the geometry appropriate to the continuously heated thruster. 

(If the cylinder is short and fat, as for an LSC wave, then radiation transport 
must also be included.) 


We can conclude that for continuously heated thrusters of oper- 
ational size (say 10 cm radius) , high partial pressures of cesium will give 
considerable amounts of radiation, which must be taken into account. For 
laboratory-size thrusters (say 1 cm radius) this radiation will be less of 
a problem, though still strong in some spectral regions. 
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To see the effect of pressure and path length, we have used 
Eq. (3-15) to find the emissivity normal to the center of the base of a hemi- 
sphere of radius I at different cesium pressures. The results are shown in 
Fig. 3-11, which corresponds for cesium to Figs. 3-1 through 3-4 for molecular 
radiation. For 1 atm pressure, the emissivities are in the 0.1 to 1 range 
for the longer lengths above 4000K, while for 0.1 atm, they are about a factor 
of 10 smaller. The highest value of pil is only about 0.3 atm ft. on Fig. 

3-10, and comparison with Figs. 3-1 through 3-4 shows how much stronger cesiiom 
radiates than the molecules considered there. The molecules are band radiators, 
with a few important bands at most, while the electrons produced by cesium 
lead to continuum radiation. This accounts for the greater importance of 
cesium as a radiator compared with the molecules considered here. 
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Fig. 3-11. Emissivity of Cesium 
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rV. THERMODYNAMIC AND TRANSPORT PROPERTIES OF 
EQUILIBRIUM SEEDED HYDROGEN* 

The LSC wave model needs the thermodynamic and transport properties 
of seeded hydrogen to high temperatures (up to 20000 K or -more). The standard 
equilibrium programs are based on the JANAF Tables and are limited to 
6000 K. In the previous work on LSC waves in hydrogen (Ref. 4) a 
rather simple analytical model of the thermodynamic properties of equil- 
ibrium hydrogen was constructed, using molecules, atoms, first ions and 
electrons. It included vibration-rotation coupling of the molecules, 
lowering of the ionization potential, and a calculated cut-off for the 
electronic partition function of the atom. The results were shown in 
Ref. 4 to be in very good agreement with the elaborate calculations of 
Ref. 32. For transport properties, the results of Yos in Ref. 33 were 
utilized. 

For seeded hydrogen the model of Ref. 4 has been extended to 
include a small amount of easily-ionizable alkali metal seed. In addi- 
tion, the trace species H was added because of its importance as a 
radiator when a low temperature donor of electrons is present. This 
extended model provides the thermodynamic properties of seeded hydrogen. 

For the transport properties, the approach of Yos (Ref. 33) was followed 
with the addition of the seed species, and an update of the cross-section 
used. 

In this section we shall present the models for the thermodynamic 
and transport properties and some results for a cesium seed, relying for 
background on Appendix A of Ref. 4. We shall see that the model of 


* We have included the H species in this description even though its 
inclusion was accomplished after the period covered by this report. 
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cesium-seeded hydrogen is quite simple to use on a computer and gives 
quite good results, as far as they can be compared with other calcula- 
tions. The thermal conductivity of cesium-seeded hydrogen does not differ 
noticeably from pure hydrogen, so that the difference in the values used 
here and those from Ref. 4-3 used in Ref. 4-1 are due only to improved 
values of cross-sections of the hydrogen species. 

4 . 1 Equilibrium Composition 

The reactions to be considered are the dissociation of hydrogen 
molecules, the ionization of the seed atom (denoted by S) , the ionization 
of the hydrogen atom, and the detachment of electrones from H . This 
comprises four reactions among seven species, as follows: 


»2 


2H 

"h''’'h2 “ 


(4-la) 

s 


+ 

S + E 

Vsp/"s ' 


(4-lb) 

H 


+ E 



(4-lc) 



H + E 

"hV"hm * 


(4-ld) 


We have written the corresponding law of mass action to the right of each 
reaction, with n^ the number densities and K the equilibrium constant for 
each reaction, which will be defined later. The subscripts H2, H, HP, E, 
S, SP, HM denote the seven species: hydrogen molecule, atom, first ion, 

electron, seed atom, seed first ion and hydrogen negative ion. There 
are seven variables in these four equations, for fixed pressure and temp- 
erature . 

There are three additional relations which make the solution for 
the n^ a determinate problem. They are the conservation of charge 
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(4-2a) 


”e ” ”hp "sp " ”hm 


the specified ratio of seed to hydrogen atoms 


' = <”s * "h ■" V" "hm' 


(4-2b) 


and the pressure 


p = nkT, n = * "s * V" "hm * "e 


= n„_ + n„ + n_ + 2n,,_ + 2n__ 
H2 H S HP SP 


(4-2c) 


For given p, T, f these 7 relations (4-la,b,c,d) and (4-2a,b,c) determine 
the solution for the seven n. . 

The complete solution is quite complex, and we will follow the 

approach used successfully for pure hydrogen in Ref. 4 by solving 

* * 

simpler problems in each of two regions T < T and T > T where in 
degrees K, 


T* = 2000 log^Q P (atm) + 7000 (4-3) 

In the low temperature region H and H are the main species and the rest 

^ . + 

are trace species, while in the high temperature region H, H and E are 
the major species. In both regions the seed is a trace, so we are re- 
stricted to seed fractions of a few precent or less, which is adequate 
for providing laser energy absorption in laser-heated thrusters. 

* 

Low Temperature Composition; T < T 

Here the main species are and H, so that reaction (4-la) is 
the' principal one. We introduce the number of pairs of H and the 
associated pressure by 
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^’o = ^H2 V^' ^ == "■ 


Then Eq. C^-la) can be written as 


( k) " k t ‘ ) 


s » 


PH 


*(%y 


which provides a solution for n as 

H 


H 

2 n' 

o 


1 + 


JEi 

V’’ 


n -h 


(4-5) 


An equation for is obtained from Eq* (4-2b) by substituting n^ 


from (4-4) , to get 


"s ' * <V * "hm> ^ - "sp 


(4-6) 


When this is put into Eq. (4-lb) , together with Eq. (4-2a) for n , a 

E 

quadratic equation for results, whose solution is 


"sP " ‘2 '•'s "hp - "hm> 


* "s" »"o + "bp + "hm’ 


1 . . v 2 

-7 (K + n - n ) 

4 S HP HM 


(4-7) 


Similar use of Eq. (4-2a) in Eqs. (41c and d) yields from the 


former a quadratic for whose solutxon is 


V ' ■ J '"sp - "hm> * 


4 '"sP - "hM>^ + 




(4-8) 
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and from the latter 


"hm * '"sp V '‘S™ * "h> 

Finally, n is related to n by n' from Eq. (4-4) as 

H2 H O 

since n^p and n^^^^ are small in this low temperature region, an 

“I 

iterative procedure is used to solve for the n^ in the following manner: 

1) Guess p ' from the given pressure p by ignoring all dissocia- 

H 

tion and ionization in Eqs. (4-2b,c) so that 
PA = P/(l + 2f) 

2) Find n /2n' from Eq. (4-5), n' fromp' in Eq. (4-4), and 

Ho O H 

n and n , the latter from Eq. (4-10). 

3) Assume n and n negligible in Eq. (4-7) and find n . 

HP HM \ SP 

4) Assume n„^ negligible in Eq. (4-9) and find n„. . 

HP HM 

5) Find n from Eq. (4-8) . 

HP 

6) Recalculate n from Eq. (4-7) . 

7) Find n from Eq. (4-2b) . 

8) Find p^ from Eq. (4-2c) and compare with given p . 

9) If P„, is not close enough top, re-guess P' as p' P/p and 

T H ri T 

repeat steps (1) - (8) until p^ agrees satisfactorily with P . 

This process converges in a few iterations to a set of with 

n given by Eq. (4-2a) . 

E 
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* 

High Temperature Composition; T > T 

Here the main species are H, H^ and E so that reaction (4-lc) is 
the principal one* We introduce the number of H nuclei and the associated 
pressure by 


2n" 

o 


H 


+ n i P" 
HP' ^H 


kT (n„ + 2n^) = 2kTn” (1 + n^/2n ”) (4-11) 
H HP O HP O 


By using Eq. (4-2a) to eliminate n from Eq. (4-lc) , replacing n 

E H 

by n", and introducing p" from Eq. (4-11), we get 
O H 



(4-12) 


A similar relation for n^^ is obtained from Eq. (4-lb) . We first 

use Eq. (4-2a) to eliminate n and Eq. (4-2b) to eliminate n in favor of f 

E S 

and n”. Then we introduce 
o 


p" = 2kTn ” (f + n /2n”) 
S O SP O 


and Eq. (4-lb) can finally be written 


( 




- K f 
s 


( 2 n ")2 

o 


HM 


) 


‘"hp ' "hm> 


2n" 

o 


= 0 


2 

K kTf 
s 



(4-13) 


(4-14) 


We now have two equations (4-12) and (4-14) for n^^ and The 

equation for n has already been expressed in Eq. (4-9) from reaction 
HM 

(4-ld) ; from Eq. (4-la) with the use of Eqs. (4-11) we find 
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^ (1 - n^p/^n”) P” ^ 

2n" 1 + n,„/2n" K kT ' 

o HP o D 


2n" 

o 


n 




HP 

2n" 

o 


(4-15) 


The pressure, from Eq. (4- 2b, c) , can be expressed as 


P = p; + + kT (1 + 2£) + kT£n^ 


(4-16) 


To obtain the n. from this set of relations we must first solve 
1 

the simultaneous equations (4-12) and (4-14) for n /2n". They are of 

HP O 

the form 


aiXi + x^x^ - b^ = 0 


^ 2^2 ^ 1^2 ■ ^2 = ° 


These can be put into the form of a quadratic for : 


J( \ - ^ V 
^\2^2^ / 


^2^1 

^^2 


(4-17a) 


which yields the solutions 


/ 

^21^ + l/(x^/x2) 


= 




(4-17b) 


With these, we can now set up a procedure for finding the n. 


1) Guess p" from given pressure p by ignoring all dissociation 
H 


and ionization in Eqs. (4-2b,c) so that 


p'^ = p/(l + 2f) 
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2) Guess p" = p - P" 

S ^ 

3) Solve for ~ ^2 Ecp . (4-17) 

with n ^ = n = 0, so that 
H2 HM 

ai = 1 + K^kT/p^ , = K^kT/p|^ 

a2 = 1 + KgkT/p^ , = KgkTf^/p^ 


4) 

Find 

/2n" and 
o 


from Eq. (4 

-15) 


5) 

Find n" 

from Eq. 

(4-11) 





o 






6) 

Find njj 2 

' ^HP' ^SP 

using n" 
o 



7) 

Find n „ from Eq 
HM 

. (4-9) 



8) 

Repeat steps (3) 

- (7) 

with the same a^ - 

and a^f but 


K 

: kT 

n n 




K — _ 

I . 

HM 

HP 


> 


^1 ~ 

p" 

2n" 2n" ' 






o 

o 





K^kTf^ 


/ 2n n 

\ 

n n 


b = 

S 

+ s 

[ H2 HM 

\ . 

SP HM 



Pc 

2n" 

V 2n" 2r." 

1 

2n" 2n" 



s 

o 

\ o o 

/ 

o o 

9) 

Find n 

S 

from Eq. 

(4-2b) 




10) 

Find p^ 
T 

from Eq. 

(4-2c) 

and compare 

with 

given p 

11) 

If p^ is 

; not close enough to p, re- 

guess 



and repeat steps 

(2) - 

(10) until p 

^ agrees satisfactorily 


with p- 







This process converges in a few iterations to a set of n^, with 

n given by Eq. (4- 2a) . 

E 
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4.2 Thermodynamic Properties 


The properties of hydrogen were defined in Ref. 4 using the same 
low temperature - high temperature approximation used here. To find the 
properties of seeded hydrogen, those properties can be modified rather 
easily to include the seed. 

The density, for example, is defined as 


p =X)min. 

where m. is the mass of the ith species. The part of the sum over the 

tj 

hydrogenic species was calculated in Ref. 4, and may be called p . To 
this we must add the contribution of the seed species n and n (ignor- 

O ^Ir 

ing the small mass of the electrons) , which is 


<"s "sp’^ ' <"S + "sp’ * "h * V* "hm’ 


Thus we see that 


p = p^ + p^ = p^ (1 + nig/m^) = p^ (1 + f Mg) 


(4-18) 


where M is the molecular weight of the seed, and f is the seed to H ratio, 
as defined in Eq. (4-2b) . 

The enthalpy of the mixture is also a sum over species where h^ is 
the enthalpy of a given species and = ’^i^i mass density of that 

species . Then 

ph 

hydrogen seed 
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We may write this, using Eq. (4-18) , as 


h 



h®) 


h" + h" 


1 + 




(4-20a) 


where the first term is the hydrogen contribution and the second is the 
seed contribution 


hydrogen 


/ H 

/p / 


seed 


H 


(4-20b) 


The pure hydrogen contribution has been given in Ref. 4, Appendix 
A. With the present notation it is 


h« = 


W [(0.7 -H 0.2^^)^ + 

H L o 


(l * 0.2 ^ + 0.4 ^ (l - 


(2 4 0.2 ^4 0.4 V);^] 


(4-21) 


Here, the dissociation temperature, 0^^, lowered ionization temperature 

0__ , vibrational enthalpy function i) , original number of hydrogen mole- 
IL V 

cules n , and hydrogen atom electronic partition function sums S and S, 


are 


,0 

^d"’h2 


.0 

*’il"h 


IL 


% = 


_ 6v/T 




®v = ^ = \2 ^ 


(4-22) 


S. = ^ j^exp (e^^/T^j^), S 2 = ^ exp (e^^/T^j^) 


D=1 


3=1 
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The species H was not included in Ref. 4. Although it is included in 
the composition here because of its importance in radiation, it is a 
minor species and its contribution to the enthalpy and specific heat will 
be ignored. Such an approximation was found excellent for pure hydrogen 
in Ref. 4. 


We can construct the seed contribution to the enthalpy in a simi- 
lar manner. For each species there is a translational contribution 
5kT/2m^. For atoms and ions there can be an additional electronic contri- 
bution ^ and for ions, their ionization energy per unit mass The 

electronic contribution is expressed in terms of an electronic partition 

function ^ as 

e £ 



kT^ 

m 


3T 


(4-23a) 


where Q is a sum over states with energy e and weight g • : 

^ V- -e.AT 

^eZ 2 ^ 9j e 3 » g,p Q kT^ 


L 

3=0 


‘e£ 


so that 


e£ 


E e.g, e 
^ 0 3 
- 3=0 


-e^/kT 


m Q 


(4-23b) 


eX, 


The enthalpies then are 

hg = 5kT/2m^ 

hg = 5kT/2mg + h^^ 

hgp = 5kT/2mg + 


(4-24a) 


(4-2413) 

(4-24c) 
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when put into the second sum in Eq. (4-20b) the masses cancel, and we have 


P-^- - 2n ) + n m h + n^ m^ (h h^^ 

1^11 2 S SP SS eZS SP S ei2,SP IS 

seed 


The factor 2 multiplying n^^ includes the contribution of the electrons 

from the seed. Note that for many seeds the electronic contribution to 

the ion enthalpy is negligible because the excited state energies e . are 

H ^ 

very large. The density of the hydrogen is p = 2n^m^, with n^ defined 
in Eq. (4-22) , so Eq. (4-20b) yields 


(4-25) 


h^ 


5kT 


2n 


2n 


SP 

2n 




Veils 

kT 


+ 0.4 


[|sp "s’''etsp^‘'is> 

2n kT 


The specific heat of the mixture is obtained from 
= 1 

H S 

obtained by differentiating Eq. (4-20a) , where h and h are given in 
Eqs. (4-21) and (4-25). 

Actually, the low-temperature, high- temperature approximation per- 
mits some simplification because for the hydrogen only the species H and H, 

+ 

or H and H appear, so that only two of the three terms in Eq. (4-21) are 

needed at any one time. The same is true of the derivatives of the n^, which 

are needed for c . 

P 
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Low Temperature Properties ; T < tT 

Here the species is ignored (n = 0) so the hydrogen enthalpy 

HP 

from Eq. (4-21) becomes 


5kT 

h = 




(0.7 + 0,2iJ,- ) 

V n 


H2 




(4-27) 


The term 1 - S^/S^ is ignored because the first terms in the sums are 
equal , and they dominate at low temperatures . There is no corresponding 

s 

simplification for h , except that for most seeds the term is 

negligible in Eq. (4-25) . 

The species ratios have already been given in the previous 

subsection, where n was called n' for the low temperature approximation. 

o o _ 

For c , the derivatives of ^ 2 :e needed; without H they can be given 

by a simplified version of the procedure outlined above for 

follows. 

The two reactions occurring here are (4-la) and (4-lb) . With 
only electrons from the seed, not from the hydrogen, these are 


H 


' Vh2 


SP 


s s 


(4-28) 


The original number of molecules, and the seed fraction from Eq. (4-2b) 
are 


n 

o 


n, 


H2 


+ n^/2 


f = 


(n_ + n__)/2n 
S SP o 


(4-29) 


while Eq. (4-2c) shows the pressure to be 

p - kT + ■>„ + "s + 2"spl - V^'’o * ^ 


(4-30) 
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We can solve Eq. (4-29) for n and n as 

H2 S 


^H2 ^ ^ ^ 

n 2n ' 2n 

o o o 


n^ 


= f - 


SP 

2n 


(4-31) 


and put them into Eq. (4-28) to get 


2n I 4n \ 2n / \ 2n 

O' o \ 0/^0 


/ n 

-i (f ._§£ 
2n V 2n 
o \ o 


By using Eq. (4-30) to eliminate the first n^ on the right we find 


2n / K kT 
o ' L D 


= 1 + 21 - 


H 

2n 


f + 


SP 

2n^ 




We can now differentiate with respect to T at constant p, re- 
membering that K and K are functions of T (to be given later) . The 

LJ 

result is two linear simultaneous equations for the derivatives of 
nj^/2n^ and which are needed for c^. These equations are 



(4-32a) 

. dlnK^T 

4p D 

Kj^kT dT 
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The solution of these two equations gives the necessary derivatives. The 

other derivatives of n„/n and n_/2n are found from Eq. (4-31) as 

Ho So 



(4-33) 


Using Eqs. (4-31) and (4-33), the 
be written from Eqs. (4-27) and (4-25) as 


H S 

derivatives of h and h may 
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Here we have defined 


0^^ = ni_h^/k and (|) = d\\) T/dT = ^6v/T 

IS SI W 


(4-35a) 


and the derivative of h « is found from Eq. (4-23b) to be 

e X/ 


dT k 


-1 


^2 ^2 
Q, 


-el 




e 


-e . /kT 
3 


(4-35b) 


When the solutions of Eq. (4-3a) are inserted into Eq. 

expression for c follows immediately from Eq. (4-26) , and h is 
P 

putting Eqs. (4-27) and (4-25) into Eq. (4-20a) - 


(4-34), the 
found by 


•k 

High-Temperature Properties: T > T 


A similar procedure is used here, where there are no H^ species. The 
hydrogen enthalpy, Eq. (4-21) , is 



(4-36) 


The reactions occurring are now (4- lb) and (4-lc) . The conserva- 
tion requires n^ = n^^ + ignoring the H species, so the reaction 

equations may be written 


(h + 
HP 




n. 


HP 


= K_ 


"h' ‘"hp 


+ "SP> 


SP 


S S 


(4-37) 
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In this regime 


"o ' ‘"h * ^ ' '"S ^ "sp’^2"o 


p = kT (n„ + 2 + 2x1^^) = 2 n>T (1 + n„^/2n^ + f + n^„/2u^) (4-38) 


H 


HP 


HP' 


SP' 


so that 


n, 


^ = 1 - 
2n 


n 


HP 


2n 



(4-39) 


The reactions, Eqs. (4-37), become 




( 1 - 


HP 


2n 


) 



n„_n__ 
HP SP 

4n 2 
o 



o o 


and elimination of the first n on the right in favor of the pressure^ from 

o 

Eq. (4-38) yields 
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The derivatives of these yield the two linear equations as 




The solution of these gives the derivatives needed for c^, along with the 
derivative of Eq. (4-~39) , which is 



(4-41) 
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JJ 

The derivative of h is found from Eq. (4-36), using Eqs. (4-39) 
and (4-41) . The result is 



- 51^ 


A 

+ 04 II* ( 

< \ 

(!± _ 

/p 

% 

^ ^ 2n 
o 


2n / 
o 



+ fl + 0 4^^^ —) 

m 


^HP \ 

V T J 

X 

\ 3T 

2n / 

° pJ 


(4-42) 


where we define 


L 

I 

j=l 


.-2 

3 exp 




(4-43) 


The expression for Oh /3T)p is the same as Eq. (4-34b) and need not be 
repeated 


Again the expression for c is given in Eq. (4-26) , when the solution 
for the derivatives from Eqs. (4-40) are inserted in Eqs. (4-34b) and (4-42) . 
The enthalpy is found by using Eqs. (4-25) and (4-36) in Eq. (4-20a) . 


Comparison of the results presented above with those given for pure 

hydrogen in Appendix A of Ref. 4 shows them to agree for no seed, f = 0 

after correction of some misprints in Ref. 4. (In Eq. (A-27b) of that 

reference 0^ should be 6^^, and in Eq. (A-29b) , the first should be 

(0 /T)^). In the case of no seed, of course, the simultaneous equation 
XL 

pairs Eqs. (4-32) and Eqs. (4-40) reduce to single equations for the derivatives 

of n /2n and n„ /2n . These latter ratios are called 3 and a in Ref. 4w' 

H o HP o 
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4 .3 Equilibrium Constants 


The quantities on the right side of Eq. (4-1) are the equilibrium 
constants for the reactions, which come from the law of mass action. For 
a reaction of the form 


A — ^ B 


this relates the partial pressure p^, p^, p^ to the temperature, the internal 
partition functions and the molecular masses. 


P P 
B C 





m. 


3/2 


(kT) 


5/2 


®iB ®iC -h^ m^AT 
— e A A 

®iA 


= KkT 


(4-44) 


Here h^ is Planck’s constant, and h^ is the heat per unit mass of A to be 
supplied to produce B + C. The ratio on the left can be converted to 
by merely changing the power of kT from 5/2 to 3/2, since p^ = nAT. This 
law of mass action can now be applied to each of the reactions (4-1) to find 
the appropriate K, which is the right side of the above expression divided 
by kT. 


Hydrogen Dissociation: Reaction (4-la) 


Here A = H^, B = C = H, and h^ = h° 


The internal partition function 


for H^ is the product of the contributions of states due to electrons, 
^^H 2 f Qr vibration Q^, and vibration-rotation coupling, Q^. 

Standard statistical mechanics gxves these as 




= 1 , Q = T/20 , e = 85.34, Q = (1-e 
IE2 ' ' r r v 


O - T + 28.45 5 5.04 E-2 E-2 

Qc - 1 + T + 1.785 E bT + ^ 
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The first recognizes that the electronic states of are too high to be 

included (greater than 132,000K) , so the ground-state statistical weight 

is used. The second is the fully-excited rotation form of Q , which is 

r 

satisfactory for high tenperature, since 0^ = 85.34F. The factor 2 
accounts for the symmetry of the hydrogen molecule . The third expression 
is the vibration of a harmonic oscillator, 0^ = 5978K. The coupling 

expression comes from Ref. 34 p.l63, Eq. (7.39), where it is given as 

and the constants are found from the spectroscopic 
constants for as explained in Ref. 4, Appendix A. 

The only internal partition function for H is due to electronic 
states, and can be expressed as 

-e A 

where is the sum given in Eq. (4-22) . The upper limit L of the sum S^, 

^(and of sums S and S defined in Eqs. (4-22) and (4-43) must be finite 

2 

since the terms grow like j for large j. Related to this is the fact that 

that energy h^ in 0 is not the usual ionization energy h^ for hydrogen, 

but a lowered value. They are both associated with the pressure of free 

electrons in an ionized mixture. A description of how to calculate L and 

6 from the number of electrons and 0_= 157,900K is given in Appendix A 
XL I 

of Ref. 4. 



Finally, ~ 52,000K = Sjj- So the equilibrium constant 

for the dissociation reaction is 

3/2 





3/2 


2 

_2eliL 


Q Q Q 

c V r 


-6dA 


(4-45a) 


"^2^^ " 52,000K 
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The logarithmic derivative of is needed for Eq. (4-32a) , in the low 

temperature regime. We ignore the derivatives of Q and and use 

Q and Q from Eq, (4-44) . 
r V 

The result is 


dJln KjjT 
dT 



(4-45b) 


Hydrogen Ionization; Reaction (4-lc ) 


^iE 

then 


Here A=H, B=H+, C=E. The internal partition functions are 
' ®iHP ~ ^ Qg2 ,H given. The equilibrium constant is 



(4-46a) 


6^ = h° m„A = 157,900 K 

X X H. 


The logarithmic derivative of this is used in Eq. (4-40a) , and is 


diln K^T 
dT 



(4 -46b) 
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Hydrogen Detachment ; Reaction (4-ld) 


Here A=H, B=H, C=E. The only new internal partition 
function is that for E~ , so the equilibrium constant is 


27Tm„k 

E 




0. = h ra A 

A AH 


= 8759K 


Seed Ionization; Reaction (4-lb) 


Here A=S, B=S+, C=E. The internal partition functions for S 
and S+ depend on the particular seed of interest. In any case 


2TTm„k 

E 


^3/2 2Qg£gp _0__/T 
0 lb 


6__ = m^/k 

IS IS s' 


The logarithmic derivative of K appears in Eqs . (4-32b) and (4-40b) . 

b 

The derivatives of Q^^can be related to ^g^^jy Eq. (4-23a) , so we may write 


d)ln KgT 


^IS ^ "^s^eilSP _ 

T kT kT 
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The seed we have used so far is cesium. For this seed the 


properties are 


®IS 45200K, 


eilSP 


= 1 


(4-49) 


eils 


S 

j=0 




-e . /kT 
e 3 


where the weights are energies can be found in the standard tabulations, 

such as the JANAF tables. We have used 23 terms in the sum, and the values 

of g.r G. and e.A are tabulated in Table 4-1. 

3 3 3 

We have now specified all the properties necessary to find the 
equilibrium composition and thermodynamic properties of the hydrogen-seed 
mixture. For pure hydrogen, they reduce to those of Ref. 4, which were 
shown there in Table A1 to agree very well with the elaborate calculations 
of Patch, Ref. 32, up to 100 atm and 20000K. 
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TABLE 4-1 


Statistical Weights and Energies for Cesiu m 



g. 

E . (cm 

e.A (k) 

i 

11 

3 

3 

0 

2 

0 

0 

1 

2 

11178 

16083 

2 

4 

11732 

16880 

3 

4 

14499 

20861 

4 

6 

14597 

21002 

5 

2 

18536 

■ 26670 

6 

2 

21766 

31317 

7 

4 

21947 

31577 

8 

4 

22589 

32501 

9 

6 

22632 

32563 

10 

2 

24317 

34987 

11 

8 

24472 

35210 

12 

6 

24472 

35210 

13 

6 

25764 

37069 

14 

26 

26616 

38295 

15 

34 

27364 

39371 

16 

72 

28446 

40928 

17 

38 

29429 

42342 

18 

28 

30172 

43411 

19 

66 

29148 

41938 

20 

96 

29679 

42702 

21 

130 

30042 

43224 

22 

168 

30302 

43598 
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4.4 Transport Properties 


The calculation of the transport properties of gas mixtures 
has been studied extensively. An exact theory is extremely complicated, 
and has been worked for only a few cases, such as ionizing argon. Most 
engineering calculations are based on various kinds of mixture rules . 

Perhaps the most successful such calculations were done by Yos (Ref. 33) , 
who presented results of calculations for air, oxygen, nitrogen and 
hydrogen; each gas was considered a mixture of its equilibrium constituents. 

Yos* work has been widely used and revised. A recent revision (Ref. 35) 
did not include hydrogen, but it did up-date some of the cross-sections 
which are needed for a hydrogen mixture. 

There is little information available on alkali metal vapor, 
especially at high temperature, and none that could be found on seeded hydrogen. 
Since the seed is a small fraction of the mixture, one might feel that the 
transport properties of pure hydrogen would be satisfactory. However, the 
seed produces electrons at low temperatures, and they are a good conductor 
of heat because of the ease with which they transport energy. So it was 
considered necessary to check the possible effect of seeing on hydrogen 
thermal conductivity. This was accomplished by using the Yos method to 
calculate the transport properties of cesium-seeded hydrogen, and comparing 
the results with pure hydrogen. This new calculation also afforded the 
opportunity to use the more recent cross-sections of Ref. 35 to up-date the 
properties for pure hydrogen. This results in significant differences 
from Yos* thermal conductivity at high temperatures. 

We will see that the thermal conductivity of hydrogen is not 
affected noticeably by seeding with a percent of cesium, especially when 
the degree of uncertainty in the cross-sections used is taken into account. 

It is therefore sufficiently accurate for the LSC wave model in seeded 
hydrogen to use pure hydrogen transport properties - 
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Thential Conductivity 


Yos expresses the thermal conductivity as the sum of three 
terms giving the contributions of the translational motion the internal 

degrees of freedom, and the change in composition due to reactions 

Each contribution is expressed in terms of the number density of the 
ith species n^, the stochiometric coefficients of the ith species in the 
nth reaction a . and a number of other parameters, as follows: 


a. . 
ID 


1 + (1 - m^/mj (0.45 - 2.54 

(1 + m . /m . ) ^ 

1 D 


(4-50) 


V. . = 

3-D 


1 

”2m.m. *1 

-h-2 

Tr(m. + m.) kT I 

Lid J 


(4-51) 


3-D- 


( 1 , 1 ) 


( 2 , 2 ) 


— V. . 7T fii-i 
3 13 


. (2) 16 V. . TT . 

' \j = 5- ‘3 


(4-52) 


The 7T Q are collision integrals, defined as weighted averages of the 
collision cross-sections: 


_ (t,s) 

TT . 

ID 


dY e Y ^ ^ (1- " cos^ X) 4ir a, ^sinX 


(s + 1) ' Tl - ^ ^ 

IS + i; . ^1 2 (t + 1 ) J 


3-D 


(4-53) 


Here (X» g) is the differential scattering cross-section for the pair 

of species i - j , X is the scattering angle in the center of mass system 
and Y is the relative velocity of the colliding particles normalized by V. . . 
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The denominator is the rigid sphere value of the double integral in 
the numerator, which is obtained by putting 47T a. . =1. 


In addition, the internal contribution involves the 

dimensionless internal specific heat, c . . and the reaction contribution 

rnt, 1 

involves the dimensionless heat of reaction A h /T, In these terms the 

n 


thermal conductivity is 


X = X + A + X 

C tr int r 


(4-54) 


where 


X 


tr 


15k 

4 



n. 

1 


a . n 
13 


^( 2 ) 

j 


(4-55a) 


X 

int 


= k 



n. c, 

1 int,i 



X 

r 


= k 



n 1 


(Ah^/T)^ 

(a ,/n.) E (a .n.-a n. ) 
nx X , nx j n] X xj 


(4-55b) 


(4.55c) 
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The sums on i and j go over all species, while the sum on n goes over 
all reactions. 


In order to evaluate this expression we need the nvimber 

densities, whose calculation has been described above, the internal 

specific heats, which can be obtained from the internal partition 

functions already given, the stochiometric coefficients, which follow 

— ( s t ) 

from the reactions, and the collision integrals tt j ^ for s = t = 1 
and s = t = 2. In these calculations we will ignore the minor species H-, 

The reactions (4-la) and (4-lb) proceed in the low temperature 
region. The Yos formulation requires them to be written as a sum on the 
left equals zero on right. When they are so expressed we see that 

^ ^1H2 = “2S ■ 1' ^2SP =-!' ^2e“ 

In the high temperature region reactions (4-lb) and (4-lc) proceed so 
that 


T > T*: a^g = 1. a^gp = -1, = -1, = 1, a3^^p = -1, a3^ = -1 (4-56b) 


The other values of a . are zero in the two regions. 

ni 


The dimensionless heats of reaction are defined as 


Ahn/T = y a . h. m.A 
ni 1 1 

X 

For reaction (4-la) this is 


Ahi 


(h.. -h ) 


kT H2 H 
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The two enthalpies can be inferred from Eq. (4-21) , or found in 
Eqs. (A-24) and (A-23) of Ref. 4. The molecular enthalpy is the 
sum of translational, rotational and vibrational contributions 


7kT 


H2 


2m 


H2 


^H2 ^ 


(4-57) 


where ij; is defined in Eq. (4-22) . The atom enthalpy is the sum of 

V 

translational , dissociation and electronic contributions 


h 


H 


5kT 

2m„ 


'IL 



(4-58) 


where S and S are given in Eq. (4-22) . When these are out into the 
1 2 ^ 

heat of reaction for reaction (4-la) we ignore the electronic contribution, 
since this reaction is only considered in the low temperature regime. 

We then find 



For reaction (4-lb) we have 


Ah^ m m 

T k ^SP k ^ 


(4-59) 


The enthalpies have already been given in Eqs. (4-24) and the heat of 
reaction becomes 


Ah 


2 


T 


1 

Oe£S 


I 

j = 0 


fill 

kT 


-e./kT 



5 

2 


(4-60) 
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For the third reaction, (4-lc) 



r % - V - 4-^ 


The ion enthalpy has translational, dissociation and ionization 
contributions 


‘HP 


5kT 


+ h: 


IL 


(4-61) 


so that use of this and Eqs. (4-58) and (4-24a) yield 



(4-62) 


These three expressions, Eqs. (4-59), (4-60) and (4-62), provide the 

heat of reaction terms for A . 

r 


The internal specific heats can be found from the 

enthalpies by differentiation. The electrons and the H+ have no internal 
degrees of freedom and so ~ ^ them. For the other species, we 

ignore the 5k/2m term in the enthalpy and form = mOh/8T)p/k from 

the other terms . 


For use of Eq. (4.53) yields 


c, ^ = 1 + $ 

int, H2 V 


where $ has been defined in Eq. (4- 35a) . 

V 


(4-63) 
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For H only the last term m Eq. (4-58) contributes, and 
differentiation yields 

/ \ 

^2 


°int, H 



(4-64) 


S, 


where is defined in Eq. (4-43) . 


For both S and S+, the only term which contributes can be 
seen from Eqs. (4-24b, c) to be h . so 


m_ dh 


eZS 


"int, S 


dT 


c . 


int, SP 


^ dh 
k 


e£SP 


dT 


(4-65) 


where the expression for the derivative is given in Eq. (4-35b) . For most 
seeds, including the cesium case calculated later, the excited states of 
S+ are so high that they may be ignored, and SP ~ 

Equations (4-63) - (4-65) provide all the non-zero expressions 

for c. ^ that are needed in A. 

int int 


Finally, there are the collision integrals tt 

- (2 2 ) ^ 
and 7T n. ' . We are considering 6 species: H,., H, H+, E, S, SP. 

Thus the matrix of the integral is 6 x 6, but because of symmetry there 

are only 21 different values in each of the two matrices. 
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The interactions between charged particles are associated 
with Coulomb forces. In Ref. 35 the Yos collision integrals for these 
interactions are updated to (cgs units) 


_( 1 , 1 ) _{ 2 , 2 ) 
TT Q. . = 7T jj. , 


Q 


c 



In MRS units this becomes 


(4-66a) 


4.38E-10 


1 

2 


An 



.82E-13 

n 


E 



0.577 


(4-66b) 


Reference 35 further suggests that for collisions involving electrons we 

use 0.6 Q , while using Q for all other Coulomb collisions. This then 
c c 

provides both collision integrals for the collision pairs: H+ - E, S+ - E, 

H+ - H+, S+ - S+, E-E, H+ - S+. 


The values for the H - E system, as given by Yos, but in 
MKS units, are 
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The other collision integrals for hydrogen species with 
other hydrogen species or electrons are also taken from Fig. 1 of Yos. 
These are H-H, ^2*^2 ^2 " latter case, the two 

collision integrals are taken to be equal. The collision integrals for 
- H+ are ignored since in our model of the thermodyraanics H 2 appears 
only in the low temperature regime and H+ only in the high temperature 
regime. 
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We have so far covered 13 interactions, leaving only the 
8 interactions between the seed species and themselves, or the seed species 
and the hydrogen species where only one of the two interacting species 
is ionized. These are: - S, H - S, H+ - S, - S+, H - S+, S ~ S+, 

E - S, S-S. There is little if any information on these interactions. The 
last three are pertinent to the pure seed, but even information on them is 
lacking for alkali metals at high temperatures. Fortunately, the seed will 
be at most, a few percent of the mixture, so rough approximations will be 
satisfactory. We took the seed atom and ion to be equivalent to the hydrogen 
atom and ion, so that we took the collision integrals for the following 
interactions to be the same. 

H 2 - S = - H; H S = H - H; E - S = E - H; 

- S " S = H - H; S - S+ = H - H+ 

The remaining three, - S+, H - S+, H+ - S, are harder to estimate. 

The latter two cannot be approximated by H - H+ because that is a charge 

exchange process with a large cross-section, while replacing one partner 

by a seed does not yield a charge exchange process. We did calculations 

—18 —20 2 

with several values, of the order of 10 to 10 m . 

Table 4-2 summarizes the source of the collision integrals. 

The curves of Fig. 1 of Ref. 33 used were read at T = 1000 (500) 10,000 
(1000) 16,000 (2000) 30,000 and interpolated linearly in log 7T fi. 

All the information is now available to calculate the thermal 
conductivity First calculations were made for hydrogen to compare with 
Yos to check the program. Calculations were made using both Yos* Coulomb 
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TABLE 4-2 


Source o f Collision Integrals 




H+ 




Fig.l Fig.l — 


E 


S+ 


Fig.l 




Est . 


H 


Fig.l 


Eq. (4-68) 
Fig. 1 


Eqs.(4-67) H-H 


Est. 


H+ 


Eqs.(4 66) 


0 . 6 * 

Eqs . (4-66) 


0 . 6 * 

Eqs- (4-66) 


E 


0 . 6 * 

Eqs. (4-66) 


0 . 6 * 

Eqs. (4-66) 


S 


H-H 


H-H+ 


S+ 


Eqs . (4“66) 


Note: Fig. 1 refers to Fig. 1 of Ref. 33. 


cross-section, which is 



and the one suggested by Aerotherm in Ref. 35, which is given in Eqs . 
(4-66) . The results are shown in Fig. 4-1 as plots of A. vs T for 
p = 1, 10, 100 atm up to 30,000K. The solid lines use the Yos expression 
for while the dashed lines use the more recent Aerotherm expression. 

The solid lines are in excellent agreement with the tables 
in Yos ' report (Ref 33 ) . The dashed lines show the effect of the more 
recent Coulomb cross-section at the higher temperatures and pressures, 
where ionized species are very significant. The curves have two bumps, 
the lower temperature bump reflecting dissociation of and the higher 
temperature bump the ionization of H. 

Calculations were then made for cesium-seeded hydrogen with a 
ratio Cs/H = 0.01. The results show only small changes from pure hydrogen 
at the same total pressure, of the order of a few percent, over a factor 
of 10 change in the value of the collision integrals which had to be 
estimated. We conclude that small amounts of cesium seed do not have a 
significant effect on hydrogen thermal conductivity, expecially in view 
of the uncertainty in the various collision integrals. 
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study of Eqs, (4-55) make this conclusion plausible. 

Addition of new species adds terms to the denominators of each component 
of decreasing each term in the sums; but it also adds more terms to 
the sums over i. These are compensating effects, and the net effect is 
small for small amounts of seed. 

The conclusion that pure hydrogen conductivity is a good 
approximation is important for simplifying the computer models of laser- 
heated thrusters. It means that can be treated as independent of 
composition, but only dependent on T and P. Thus we can input ^ as a 
table independent of seeds, and interpolate, rather than having to change 

A for each seed fraction. 

C 

The calculations by Yos ' method described above are limited 
to T > lOOOK by the collision integrals given in Fig. 1 of Ref, 33, 

Below lOOOK there are only two species, and S, so the calculations are 
quite simple, having only A^^ and components. It can be verified that 

at lOOOK also there is only a few percent difference between pure and 
Cs/H = 0.01. In fact, a good fit to the calculations at lOOOK and 2000K is 

/ ( 1 + 0.02f) (4-70) 

which shows the small effect of small seed fractions f. 

The pure hydrogen expression for T < lOOOK can be 
using the value A^ = 0.434 W/mK calculated at lOOOK, and the 

W/mK given on page 574 of Ref. 36 as experimentally measured, 
fit gives 

= 3.67E-4 T + 6.69E-2 (4-71) 

for pure hydrogen, T < lOOOK. 


obtained by 
value 0.177 
A linear 
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Viscosity 


Yos also gives an expression for viscosity in Ref, 33. 
in the same terms as his conductivity expression. It is 


y =X;in.n./^ 

i j 


^( 2 ) 
ir j 


(4-72) 


(2) 

where A. . 

1/1 

Eq. (4-52). 


has already been defined in terms of a collision integral in 


We can expect a smaller alkali seed fraction to have a larger 
effect on y than on because of the large mass of the seed atom compared 
to hydrogen. (The molecular weight of cesium for example is 133.) Cal^^ 
culations bear this out. At low temperature, up to 4000K, the viscosity 
of Cs/H = 0.01 is as much as a factor of 2 larger than that of pure 
This difference is reduced to only a few percent in the middle range of 
temperature but begins to become significant again at temperatures near 
20,000K. These effects will have to be considered if viscosity becomes 
an important parameter. 



V. DEVELOPMENT OF THRUST CHAMBER DESIGN COMPUTER PROGRAM 


The two inodes of heating the propellant being considered necessitate 
two different computer programs. In the continuous heating mode, both the 
low and moderate temperature absorbers are present. Laser energy absorption 
occurs continuously from the input temperature to the high temperature at 
which all the laser energy has been absorbed. The axial gradients are not 
large, so viscous and heat conduction effects in the axial direction can be 
neglected. For this case we need a flow model which includes radial grad- 
ients, laser absorption, equilibrivim chemistry, and radiation losses. 

In the discontinuous heating mode, there is no low temperature ab- 
sorber. Absorption occurs in an LSC wave standing in the flow channel, and 
the low temperature heating is accomplished by the LSC wave mechanism of 
forward energy transport by radiation and conduction. Here we need a flow 
model which includes axial conduction and radiation transport, laser absorp- 
tion, and radiation losses. The radial gradients are not important here, 
except possibly near the edges of the heated zone. 

In previous work performed at PSI, (Ref- 4) , a one-dimensional LSC 
wave model for pure hydrogen was constructed. That model is the basis for 
the seeded hydrogen L^C wave thruster being considered here. The additions 
needed to include seeding are in the thermodynamics, the thermal conductivity, 

f 

the laser absorption, and the radiation transport and losses. The models 
needed for the first three additions have been discussed in the previous 
sections of this report. The radiation transport and loss model for seeded 
hydrogen has not yet been developed. Since the LSC wave model for pure 
hydrogen is thoroughly described in Ref. 4, and the version for seeded 
hydrogen is not complete yet, we will not discuss the seeded LSC wave model 
in this report. 

The continuous heating flow model has formulated almost completely, 
and will be described here. The resulting computer program is based upon 
an axisymmetric wake flow program which was available at PSI. To adapt this 
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program to the continuously-heated thruster we had to add laser energy ab- 
sorption, radiation losses, and provide for the presence of outer walls. 

5. 1 Equations of Motion 

We consider flow in an axisymmetric channel, where the axial gradients 
are small compared to the radial gradients. This enables us to make the 
boundary layer approximation that the viscous terms need only include the 
radial gradients, not the axial ones. Let us denote the axial and radial 
coordinates by x,r, and the axial and radial velocity components by u,v. The 
mass conservation equation is then 


3rpu 3rpv 

9x 3r 


(5-la) 


In the usual fashion, the boundary layer approximation applied to the radial 
momentum equation tells us that the pressure is constant over a cross-section, 
so p = p(x). Then the axial momentum conservation equation is 


9u 


5u 


dx r gr 




aja 


) 


(5-lb) 


The energy conservation equation must include the radiation flux 
terms, which will account for both the laser radiation and the radiation 
from the gas. Written with enthalpy h as the dependent variable, and with 
a radiation flux vector whose axial and radial components are denoted by 
and S^, the equation is 



pv 


^ 1 
9r r 


9rS 



9s 

X 

9x 


u 


dx 


+ y 




9r 



Ep 


D h 
1 1 


9r 


Note that is the frozen thermal conductivity, not including the reaction 
contribution, which is the summation term in the bracket. Thus is the sum of 
the translational and internal contributions, as given in Eqs. _(4-55a). and (4-55b) 
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The term in square brackets on the right is the heat flux in the radial di- 
rection.' It is composed of the conduction part proportional to the tempera- 
ture gradient, and the diffusive part, which involves the enthalpy per unit 
mass h^ of the i-th species, a diffusion coefficient for that species, 
and the gradient of its mass fraction C^. The form of the diffusion term 
is the simplified one usually used, where the diffusion velocity has been 
expressed by Pick's law as 9c^/9r, to avoid the considerable complication 
of considering multicomponent diffusion. 

The heat flux term can be simplified by introducing the relation 
between the mixture enthalpy h and the species enthalpy and mass fraction. 

h = Vh. C 
^ 1-1 


If we differentiate, and recognize that h^ depends only on T, we find 

dh ='^.dC. + c dT, c ^ = Yc C. , c = dh./dT 
^1 1 pf pf -^p. 1 p. 1 


The quantity c^^ is the frozen specific heat of the mixture, which is just the 
weighted sum of the specific heats of the individual components. 

By replacing the temperature derivative by the enthalpy derivative 
in the last term of the energy equation, it becomes 


9h ^ 


3h 1 , dp / 3 uY 

P'" 9r r 9x 3x dx ^V9rj 


r 9r 


Cf 


pf 


^ +T (Le.-l)h. ^ 

9r ^ 1 1 3r _ 


(5-lc) 


where the Lewis number is defined as 


Le. = pD. c ^ / X 
1 ^ 1 pf ^ Cf 
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which is the ratio of diffusion of species to diffusion of heat. This 
equation shows the considerable simplification which occurs if Le^ = 1, 
when the species concentration gradients do not need to be included. 

The three conservation equations determine the three independent 
variables \x, v and h in the two-dimensional x, r coordinate system. This 
formulation is appropriate for a gas mixture in chemical equilibrium, when 
the chemical state is determined by the specification of two thermodynamic 
variables, such as T and p, plus the mixture composition in terms of the 
fraction of each kind of chemical nucleus contained in the mixture. It 
is also possible to formulate the problem in terms of non-equilibrium 
reacting chemistry; such a formulation was contained in the original wake 
program from which the present program was adopted. However, it was found 
that the laser energy addition caused problems with the nonequilibrium 
chemistry calculation. Since nonequilibrium was not considered to be 
required for the laser-heated thruster application, at least during the 
heating phase of the flow which is being addressed here, the program was 
converted to the chemical equilibrixm formulation given above. 

Part of the radiant energy flux vector ^ is the laser beam. We 

will take a parallel beam travelling parallel to the x axis, for which 

S - I and S =0. Then the laser contributes a term 9l/9x on the left 
x r 

of Eq* (5-lc) , but the laser energy is absorbed according to 

9l/9x = -k I (5-2) 

Xj 

where k is the laser absorption coefficient discussed in Section II. Thus 

li 

the laser provides a term k I on the right side of Eq. (5-lc) . (Note that 

li 

we may permit I to depend on both x and r, even though it is only propagating 
in the x direction.) The remaining part of the radiant energy flux vector 
is due to radiation transport in the gas and will be denoted by ^ , and 
the terms on the left of Eq. (5-lc) containing this vector will be carried 

p 

as V • S , the divergence of the vector. 

The final form of the energy equation is then 



PuH . pv|| * VS’> = u|| . p(|^)>V 


r 9r 


^ Z (Le.-l)h. ^1 

Cpf |^8r 1 1 3r J 


(5-ld) 


These cylindrical conservation equations are singular on the axis 
r = 0, and must be expressed in limiting form there. On the axis, symmetry 
requires that v = 0 and 3u/9r = 0. Near the axis a Taylor expansion then 
gives 


V = 




(5-3) 


where the subscript a means on the axis. For the variables h and re- 
lations like those for u apply. These relations enable us to express Eqs. 
(5-la, 5-lb and 5-ld) on the axis. All the terms with the operator v9/9r 
vanish. The viscous terms are the ones with the factor 1/r. For the momen- 
tum equation, the viscous term is^ 

i ^ u — = n + H ^ + 9y 9u 

r 9r ^ 9r ^2 r 9r 9r 9r 

9r 


The first two terms have the same limit on the axis because of Eq. (5-3), 

2 

while the last term is of order r there and may be ignored. So we find that 


( 


1 

r 


9r 


ry 



a 


Exactly similar treatment can be accorded the last term in Eq. (5-ld) . The 
second term in the mass conservation is treated similarly (after division by 
r) using the expression for v in Eq. (5-3) . The resulting equations on the 
axis become 
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(5-4a) 



+ 



0 




pu ^ + V ■ S . 


k I + u ^ 
L dx 


+ 2 


Cf 

c £ 
Pf 


'~ 3^h 

3r^ 


E 


3^C. 1 

.-1) h. 

" " 3r' J 


(5-4b) 


(5-4c) 


5.2 Transformation of Equations of Motion 

The numerical solution of the conservation equations is accomplished 
by using the stream function ¥ as one of the independent variables, instead 
of r. This satisfies the mass conservation equation exactly by setting 

r 


-I 


pur = |t- ^ pvj- = > p!. ^ ijf = / purdr, = 2 / d'i'/pu (5-5) 


The stream function has the significance of being the' mass flow through a 
cross-section of radius r, except for a factor of 2 tt. We now transform 
from X, r to C, 4', where ? = x, but is given a different name to avoid con- 
fusion. The transformation laws are 


3_ 

3r 


. 3 9 _ 9 9 

3'i' ' 3x 3? ■ a'l' 


(5-6a) 


V = u(3r/3x) 


'l'=const 


12 6' 



Thus the convective derivative operator becomes 




and the viscous operator becomes 


1 

r 


9 / . 9 


= PU^ 


2. 9 

r fpu ^ 


(5-6b) 


where f is any function* 

^ Using these transformation equations, we find the momentum and 
energy conservation expressions Eqs. (5-lb and 5-ld) to become 


9u _ 

9? "" pu d^ 


(pur^y If) 


(5-7a) 


9h , . 1 dp . . 2.. /9uf 

+ = + - + pur y j 


9C pu pu p 


9T 


2, 

pur 

c 

9c ~\ 

^ + ^(Le -l)h^ 

Pf 

^ — i 


(5-7b) 


The corresponding equations on the axis are to be found from Eqs . 
(5-4b and 5-4c) . We need to evaluate the second radial derivatives in terms 
of From Eq* (5-5) we see that near the axis 


Y = p u r^/2, r = (2'l'/p u 
a a a a 
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Then from Eq. (5-6a and 5-6b) we find 


^ = p u r 
9r a a 


9H' 


(2p u H') 




94' 


9 

9r 


pur 
^a a 



u r 
a a 


94'/ 


2p u /r 

a a 



= 2p u 
a a 



^ i 3 \ ^ 

2 94'y ~ Pa'^a 94' 


The last step follows because 4' = 0 on the axis . Thus the second r 
derivatives in the axis equations become first 4' derivatives • Then 
Eqs. (5-4b and 5-4c) on the axis transform to 


9u 

9C 


1 

pu 


d? 


+ 2y 


9u 

94' 


9h 

9? 


+ 


V.S 


R 


P 


1/1 

L + i ^ 
pu p d? 


+ 




9c * 

"E'l-i-l'hi# 


(5-8a) 


(5-8b) 


We now have two equations in the unknowns h, u to be solved, when 
expressions for I, V • ’'cf /Cp^ and dp/d^ are given. 
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5.3 Wall Shape 


The use of the ^ coordinate instead of r means that r has to be 

calculated by Eq. (5-5 ) p once the solution has been obtained at a given 

cross section x. Since u = 0 at the wall, r = r , the integration must be 

w 

done carefully near the wall. The behavior of pu near the wall can be 
obtained from the transformation Eq. (5-6a) , which shows that 

9pu ^ iPH. = £ 9(pu)^ 

3r ^ av 2 34' 


Since we expect the r-derivative of pu to be finite at the wall, while u 

is zero there, it is clear that the ^-derivative of pu must be infinite, 

2 

but the H^-derivative of (pu) is finite. Thus we may write a Taylor series 
(pu)^ = (4'-4'^) [3(pu)^/34']^ 

near the wall. The relation between r and from Eq. (5-5) can be written 



Using the expansion for 


(pu) on the left, we find 


2 

r 

w 


+ 4 


4* -4' 
w 


[-3(pu) ^34'], 


(5-9a) 
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where we have recognized that ^ and the gradient is negative since 

pu decreases to zero at the wall. 

A finite difference expression for this equation for r is 

w 

obtained by putting the derivative equal to the difference quotient 
from r to r^, remembering that pu = 0 at the latter location. Then 
Eq. (5-9a) becomes 

r^ = r^ + 4(H^ --'F)/pu (5-9b) 

w w 


This can be applied to find r^ at when r is known at by using pu 
at T. It carries the integral of Eq. (5-5) from 'F to 

The fact that the wall shape is determined as a result of the 

calculation, instead of as an input, means that the calculation is an 

indirect one. However, there is a parameter available to control the 

wall shape r (x) , since we have not yet specified the oressure gradient 
w 

dp/dK- In the stream function formulation, we are free to specify dp/dK 
to adiust the wall shape to a desirable one. Since the flow is heating, 
the density will tend to decrease, and the streamlines will tend to diverge, 
increasing the channel size. To counteract this trend, if we want to main- 
tain an approximately uniform channel area, we may increase the flow velocity 
u. Reference to the momentum equation (5-7a) shows that we then need to 
maintain a negative pressure gradient dp/d?. How to specify dp/d? so as to 
obtain a desired wall shape is an art that will have to be learned as ex- 
perience with the flow computer program progresses. 

5.4 Laser Beam Shape 

Absorption of energy from the laser beam is a heating term k I in 

L 

the energy equation (5-8b) . The intensity I is determined by the absorption 
equation (5-2) . Notice that I is a function of x, r and not ?, ¥ since 
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it does not follow the streamlines. Thus the solution of Eq. (5-2) is 
computed in x, r coordinates, and then the relation between and r at a 
fixed X = C is used to convert I (x,r) to 1 (C, H') for use in the energy 
equation. Conversely, k is a function of the thermodynamic variables, 

li 

which are defined in the ^ system, and it is transferred to the x, r 
system for use in Eq. (5-2) . 

An initial profile, at x = 0, is needed for I = to start 

the integration of Eq. (5-2) . There are many profiles that could be 
chosen, such as a uniform (top hat) profile or a truncated Gaussian. We 
have chosen to use the far-field pattern of a beam made by uniformly 
illuminating a circular aperture. This profile is 


I. (r) 

X 


I 

o 



(ar)/ar 


2 


(5-10) 


where is the first order Bessel function, is the intensity at r = 0 

and a is a parameter giving the width of the beam. 


to be 


The power P in this beam up to any radius r is found by integration 


P 

r 



4771 

27Tr I^(r)dr = — ^ 
a 


|l-J^ (ar) 


(ar)j 


(5-11) 


2 

The total power to r = is 47Tl^/a , which shows that 2/a may be interpreted 
as the radius over which a uniform beam of intensity would have the same 
total power as . 

The parameter a may also be related to the first zero of which 
occurs at 

ar^ = 3.83171 


Thus we can write I. as 

X 


I. 

X 


I 

o 



y = 3.83171 r/r^ 


(5-12) 
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where the parameter is now the radius at the first zero of I . 

In order to allow for the presence of a buffer gas near the wall, 
where no laser heating takes place, we also introduce a cut-off radius 
r^, beyond which the laser intensity is zero. Then the initial laser 
profile is 


I. = I 
X o 

= 0 



r < r 

c 


r > r 

c 


(5-13) 


This provides three parameters with which to specify the incoming laser 

profile. I is the center-line intensity; r_ is the location of the first 
o 1 / 

zero, which specifies the width of the profile? r gives the radius beyond 

c 

which there is no laser beam. 


in Eq. 


The laser power entering the thruster is found by putting r = r^ 
(5-11) . By defining y^ as in Eq. (5-12) with r = r^, we can write 


P = 4it I 1-J^ (y ) 
c o c [_ o c 




)/v 




3.83171 



(5-14) 


If desired, can replace one of the three parameters I^, r^, r^ in speci- 
fying the initial laser profile . 

The distribution of intensity, power and area is shown in Fig. 5-1 

as a function of r/r_ , up to the first zero r = r. . I. /I is obtained from 

1 1X0 

Eq. (5-12) . The power is obtained from Eq. (5-11), normalized by = P(r=r^), 

which is found by recognizing that J. (ar, ) = 0, and J (ar^ ) = 0.40276 so 

IX o 1 

P l-J^ (y) - jJ(y) 

r _ o 1 

P ” 0.83778 
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i Power for the Incoming Laser Beam. 
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. 2 2 
The area is just A = 7Tr and at it is A^ = TTr^ so that 

A/A^ = (r/r^)^ 

It can be seen from Fig, 5-1 that the intensity drops off rapidly 
with r, and 90% of the power is inside r/r^ = 0,65, The region outside 
this radius is nearly unheated. In a sense, this laser profile automatically 
provides a practically unheated buffer zone if the cut-off radius is taken 


5,5 Radiation Model 

The term V • ^ in the energy equation (5-8b) is the flux of radi- 
ation. What is needed for this term is an engineering model of the radi- 
ative flux. Such a model can be based on the rather detailed information 
about absorption coefficients presented in Section III. However, it is 
not feasible to use complete spectral information because integration 
over the wavenumber range for all calculated points in the flow field 
would be very time-consuming. 

Although computation of radiative transfer within a gas is an 
extremely difficult nxamerical problem, there are several features of the 
laser- heated thruster which suggest a simplified model may be applicable. 
First, the shape of the laser beam (Fig. 5-1) suggests that the temperature 
at the center of the channel will be much hotter than that in the outer 
region near the cooled walls. Therefore, the radiative losses will be 
most important in the high temperature central region; it is reasonable 
to sacrifice accuracy in the outer region in favor of a more realistic 
description near the center. 

Second, the absorption coefficients increase rapidly as the tem- 
perature increases, so they are several orders of magnitude larger in 
the hot region than in the cold region. Therefore, if re- absorption is 
unlikely to occur in the hot region (optical depth less than 0.5) it is 
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also unlikely to occur in the cold region, because the decrease in ab- 
sorption coefficient with decreasing temperature is much greater than 
the increase in size of the cold region over the hot region. On the other 
hand, if the hot region is optically thick so that re-absorption must 
be considered, the radiation trapping there reduces the radiation loss 
and ca\ises the temperature to remain rather constant. Outside the hot 
region, where the gas is optically thin and losses occur, the temperature 
is expected to drop rapidly. 

This physical picture suggests the use of a thick- thin model of 
the radiation flux, in which each region is characterized by absorption 
coefficients and a thickness, and depending on the resulting optical depth, 
is described by either an optically thick (re- absorbing) or an optically 
thin (emitting) flux model. The physical picture just described also 
suggests that the transition from one description to the other is ex- 
pected to occur .abruptly, and the results will be insensitive to the de- 
tails of the transition. 

In order to avoid dealing with the spectral variation over the 
whole range, a band model is used. The spectral information available 
from the detailed description in Section III is examined, and used to 
identify a few important bands, for which approximate absorption coef- 
ficients are constructed. Integration over these bands can then be per- 
formed analytically. 

The result of combining the thick-thin description and the band 
model is a rather simple engineering model of the radiation loss which 
is easily implemented in the flow field code, and yet still should give 
a description of the radiation flux within engineering accuracy. 

In the region where a radiation band is optically thin, emission 
dominates over absorption. This approximation leads to an expression 
for the radiation source V • S involving the Planck mean absorption coef- 
ficient, and can be expressed as (Ref. 37) 

’"to V (5-15) 
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where is the blackbody function defined in Eq. (3-5) , and the integra- 
tion is taken over the thin bands. 

In the region where a radiation band is optically thick, the log- 
arithmic rate of change of with respect to optical depth is small. 

The only readiation which reaches a given point comes from points where 
B^ does not differ much from its value at the given point. This approxi- 
mation leads to the Rosseland mean absorption coefficient, and to an ex- 
pression for the radiation flux vector as 


S 


R 



dB 




k dT 
0 ) 


do) 


where the integration is taken over the thick bands. To form the divergence 
of S we recognize that only the gradient in the r direction need be kept, 
since we have taken the axial qradients to be small. Thus 


(V 



1 

r 3r 



4TT ^ 3t \ 

3r 3r \ R ^ 3r/ 


(5-16a) 


--f - 


dB 


0 ) 


thk O) 


dT 


do) 


(5-16b) 


This is often called the radiation conduction approximation because the 

form of Eq. (5-16a) is exactly the same as the conductive and diffusive 

terms in the energy equation. In fact, A plays the role of a radiative 

R 

thermal conductivity. 

The sxom of Eqs. (5-15) and (5-16) provides the V • £ term in 
the energy equation. We need now to specify the bands, the absorption 
coefficient for each band, and some criteria for deciding when and how 
a band changes from thick to thin at different points in the flow field. 


the 


For each band we will specify the absorption coefficient k 
form 


1 


in 



-hpC0)/kT 


(5-17) 
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where we have separated out the stimulated emission factor. The will 
be independent of the wavenumber, since it is appropriate to an average over 
theband. Thus the integrations indicated in Eqs. (5-15) and (5-16) can 
be performed in general over the band limits O)^ to Using the ex- 

pression for from Eq. (3-5) , we can write Eq. (5-15) for a band as 


(V • s^)^. . = STTh c“k 

— thin P 


to 


2 -hpC(o/kT 

dto to e 


(0. 




(x + 3x + 6x + 6) 


] x 

" 


(5-18) 


where x = hpCto/kT. A similar use of Eqs. (3-5) and (5-17) in Eq. (5-16b) 
shows that for a band 


X 




= 15 0 ^ r 

5 * I 

IT k. 


u 4 -X , 

X e dx 


(1 - e‘^) 


This integral is not expressible in finite form, but can be integrated 
by parts to yield a form easily calculated. The result is 


ii 1 
1 


(1 - e ) 


- 1 


^ 3 -X 

2 X e 


1 - e 


“X 


+ 2x^ (x + 3) jin (1 - e ^) - 6 ^ 


-nx 


n=l n 


(5-19) 


2 (1 + i) (x + i) + x^ 
n n 
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Band Model Example 


In order to try out the band model, a model was constructed for 
a particular case and implemented in the computer program. The case chosen 
was a mixture H^/H^O/Cs = 0,945/0.05/0.005 at p » 30 atm, for a radius 
of about 12 cm. The detailed spectral absorption coefficients were cal- 
culated for this case by the methods described in Section III. The ab- 
sorption coefficient for the equilibriim mixture is shown in Fig. 5-2, 
as a function of co, for T = 2000, 4000 and 6000 K. The components of 
this absorption coefficient are presented similarly for the same three 
temperatures in Fig. 5-3. The components are that of hydrogen molecule 
vibration, KV, hydroxyl KOH, water three cesium contributions 

KCS and the four hydrogen atom contributions KH. Another important in- 
gredient of the model is the normalized blackbody function 

4 

RE = B /or 

(a) 


which is presented in Fig, 5-4 for the three temperatures. It provides 
the limits of the radiation at each wave number, and shows where radiation 
cannot be important regardless of the size of the absorption coefficient. 

The band model was constructed after study of these plots to ex- 
tract the important part of the radiation. It was concluded that the 
radiation for this case could be expressed to sufficient accuracy by the 
use of three bands (with O) in cm ^): 

Band 1: 0 < 0 ) < 0.8E4 

Band 2: 0.8 E4 < 0) < 2E4 (5-20) 

Band 3: 2E4 < o) < 

In the first band radiation from cesiiam and from electron- atom hydrogen 
Bremsstrahlung are important. In the second and third bands radiation 
from cesivmi and free-bound H are important. The absorption coefficients 
for the three bands are expressed as 
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Fig. 5-2b Absorption Coefficient as a Function of Wavenumber for 


H /H 0/Cs = 0.945/0.05/0.005. T = 4000 K, p = 30 atm. 
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Fig. 5-2c Absorption Coefficient as a Function of Wavenumber for 
H 2 /H 2 O/CS = 0.9'45/0. 05/0. 005. T = 6000 K, p = 30 atm. 
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Fig. 5-3b Component Absorption Coefficients as a Function of Wavenumber 
for H^/H^O/Cs = 0.945/0.05/0.005. T = 4000 K, p = 30 atm. 
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Fig. 5-3c Component Absorption Coefficients as a Funct ion of Wave number 
__ fbr H /h.,o7Cs^“='07^570,0570 T^^ GOO^ K, p = 30 atm. 
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(5-21 a) 


1.75 E-17 n +4.2 E-17 n , 
Cs HM 


(5-21b) 


k* = 1.0 E-17 n* +2.5 E-17 n„. 

3 Cs HM 


(5-21c) 


The cesiim number densities are those of certain excited states, which 
are expressed through an approximation to the cesium partition function 


(5-22 ) 


^ ^ -16616/T , „ -20950/T 

Q„=2 + 6e +10e 

’^Cs 


-32540/T -38130/T 

+ 20 e + 266 e 


In terms of and the cesium number density n^^ the excited state 

densities are 


Cs 


-38130/T 

266 e 

Cs ^Cs 


(5-23 a) 


** -32540/T 

n„ = 20 e /Q. 


Cs 


Cs' ""Cs 


(5-23 b) 
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. (. 


n! = (6 . 10 „ /Q, 


) 


Cs ""Cs 


(5-23 c) 


The cross-sec±.ion in Eq. (5-21 a) is the E-H Bremsstrahlung cross-section 
from Eq. (3-18) evaluated at (a) = 4500 cm which can be expressed as 


o ( 0 ) = 4500) = 1.08 E-43 T^^^ (1 + 3237 /T) e ^ (5-24) 

EH 
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where ^ is given in Eq. (2-6). This cx)mpletely defines the absorption 
coefficients for this band model example, which are to be used in Eqs. 
(5-18) and (5-19). 

The contributions to the absorption coefficients defined above 
come from cesium, H, H and electrons. The spectra of Fig. 5-3 also show 
a contribution from water in a limited range, with considerable structure. 
This contribution is always thin, and only important at low temperature. 
Rather than treat it as an absorption coefficient, therefore, we have 
parametrized it directly as an added contribution to the emission for 
the thin case: 


(V 



thin 


1.85 E-21 


T 


2.9 



(5-25) 


Thick-Thin Criteria 

At any point in the flow field, a given band is taken to be thick 

or thin, so the radiation term is either given by Eq. (5-16) or by Eq. 

(5-15). We need a criteria for the transition from thick to thin. This 

criteria involves the product of the absorption coefficient and some 

typical length, r . 

R 

Since the absorption coefficient is of the form of Eq. (5-17), 
we need to evaluate the stimulated emission factor at some mean band 
location. Examination of the spectral absorption coefficient curves 
suggests that for the purpose of choosing a band to be thick or thin, 
we can use 

-1 * -6475/T 

Band 1: 0) = 4500 cm , k^ = k^ (1 - e (5^ 

-1 * -20143 /T 

Band 2:0)= 14000 cm ,k^ = k^ (1 - e (5 

Band 3: u = 23875 an“^,k 2 = k* (1 - . (5 


26a) 

-26b) 

-26c) 
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Various methods for making the thick-thin transition were tried. 
In order to avoid an abrupt change from thick to thin we settled on 
a gradual transition, with an interpolation function. 


k r <0.2 Thin 
1 R 

k. r >5 Thick 
1 R 


0.2 < k^ r^ < 5: Interpolate . 


The interpolation was performed using 


(5-27) 


V • = (V • + f(k 

— — thxn 




(V • - (V • 


In 






(5-28) 


1 + (k.rj^)' 


The function f is 0.038 at the thin boundary and 0.96 at the thick boundary. 
As discussed above the details of this transition are not expected to ef- 
fect the results significantly. 


The choice of the length r was based on the radiation changing 

1 

by a factor of 2, or the temperature changing by 2^ = 0.841. The tem- 
perature profile near the channel centerline was fitted by a parabola 
centered on T (^ = 0) , and passing through the nearest calculated point 
to the centerline at which the temperature was down below 0.841 T = 0) , 
call it T*, at which r = r*. Then we chose 


r 


R 


1/2 

0.159 T (¥ = 0) (r*)^ 

T('l' = 0) - T* 


(5-29) 


which is the point at which T = 0.841 TC'l^ = 0) on the fitted parabola. 

These relations now completely define the way to evaluate the 
radiative loss term V • S . 
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5.6 Heat Flux to the Wall 


The wall heating is an important parameter in the design of 
a laser-heated thruster, and its expression merits special mention. 

It is convenient to consider the wall heat flux in two parts: first the 
convective and diffusive part, and second the radiative part. 

The convective and diffusive part has already been expressed 
in Eq. (5-ld) as the coefficient of r in the last term. 

The minus sign is introduced so that is positive when heat flows 
into the wall (negative gradients) . Use of the transformation deriva- 
tive in Eq. (5-6a) enables this to be expressed as 


■3c - - o 


C f 


pf 


— + V 


(Le. - 1) 
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8c. “ 

V, 1 

X 


(5-31) 


The numerical evaluation of q^ at thewallis not straightforward 
because u = 0 there while 3/3^ is infinxte, as already discussed in 

2 

connection with the wall shape. However, if we convert the 3/3r to 3/3r 

by multiplying Eq. (5-30) by 2r, we can then use Eq. (5-9b) to convert 
2 

3/3r to finite dxfference form, and get 


'Ic = c 
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h - h 
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(C. - C.) 

xw X 


_£U_ 


4('i' -’i') 
w 


(5-32) 


where the quantities without the w subscript are to be evaluated at the 
first ^ station away from the wall. 

The radiative heating contribution has thick and thin parts, 
according to the model described above. Reference to Eq. (5^I6a) shows 
the thick part to be 


151 


/ 



R 

‘^thk 


(5-33) 


4tt a 9t 
3 R 9r 


4tt 

3 


purX 


9t 

R 9^ * 


The same conversion to 9/9r produces a finite difference form 
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(5-34) 


The thin part of the radiative flux is expressed in terms of 
R 

V- • £ directly by Eq. (5-15). Thxs can be converted to a flux by in- 
tergration over r. 
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We have constantly ignored the axial fluxes, and will do so 
R R 

If we call (S ) = q , . f and recall the expression for V • 

r w thin 

( 5 - 1 5), we have 


here also. 
£ from Eq. 


R 1 

^thin r 

w 



(5-35) 


This integral can be done in r space to avoid the complications of 
u = 0 at the wall. 

The interpolation scheme described above requires that the in- 
terpolation function f (k. r ) be included, so the final expression for 

i R 

the radiative heating is 
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(5-36) 
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The sum of Eqs. (5-32) and (5-36) gives the local wall heat flux, 
which can be integrated on x to give the total flux up to any station x 
as 


2w 



(q^ + q^) • 


(5-37) 


5.7 Numerical Method 

The equations to be solved are basically the axial momentum equation 
(5-7a) , the energy equation (5-7b) , and the laser absorption equation 
(5-2). The equations are parabolic, since they are first order in C 
and second order in Y. Thus a marching scheme can be used. Knowing 
all conditions at a particular value of we may march forward to the 
next value. The process is begun with a given set of conditions at 
the initial station ^ = 0. 

The dependent variables to be solved for are the axial speed 
u, the enthalpy h, and the laser intensity I. The pressure distribu- 
tion is input. When p and h are known, the temperature and composi- 
tion are found from the hypothesis of chemical equilibrium. 

The non-linearity of the equations, and the coupling between 
u and h require an iterative solution. To advance a step in Eqs. 

(5-7) are solved for u and h by the finite-difference Crank-Nicholson 
implicit scheme. This requires the values of k , I and other properties 
at the new ^ station, which depend on the given p and an unknown T. 

To start the advance, a T profile is guessed, based on the profile at 
the present ^ station. This is used to calculate the needed fluid prop 
perties and I. Then the values of u and h are found by Crank-Nicholson 
at the new station. With this h and the given p, we find T and from 
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equilibriiHcv. If this T does not agree with the guessed profile, it 
is ^used to re-calculate the fluid properties, and the equations are 
again used to find new values of u and h at the new station. This 
process is repeated until the T profile used for the properties agrees 
with that found from the calculated values of u and h to a satisfac- 
tory accuracy. This then completes the calculation at the new station, 
and we then advance again by repeating the iterative process. 

In the laser-heated flow being considered here, the properties 

which have the main influence on the flow are k and I, through the 

L 

first term on the right of the energy equation (5-7b) , which is the 

laser heat addition term. Since k is a strong function of T and the 

L 

number densities n = p C. /m , it is the dominant effect on the tem- 

1 11 

perature profile, and controls the iteration process. 

The computer code which embodies this numerical method is called 

LHTE. 


5.8 Numerical Example 

; 

We have calculated a numerical example using the code LHTE 
for 10.6ym. Because the radiation model had only been developed for 
p = 30 atm and a mixture of H^/H^O/Cs = 0.945/0.05/0.005, this was 
the case chosen. The initial profile of T was uniform at 1000 K, and 
the wall was held fixed at that temperature. At this pressure and tem- 
perature the initial density is 1.27E-3 g/cm^ for this mixture, whose 
molecular weight is 3.47. The centerline laser intensity was taken 
as lo"^ W/cm^ = lO^^erg/cm^s and the power as 10^ W. For the beam shape 
we took the cut-off radius to be the same as r^, the radius to the 

first zero. Then Eq. (5-14) determines r = 11.81 cm. We chose a 

^ c 

case which would have, if there were no losses, an I = 2000, so 

sp 

the kinetic energy in the exhaust would be 1.92E-12 erg/g. Since this 

6 13 

is also the power per unit mass flow, a power of 10 W (10 erg/s) 

leads to a mass flow rate of 5.21 g/s. We took a uniform velocity 

profile, and set the initial wall radius r to r =11.81 cm, so the 

2 w c 

initial area was 438.2 cm . Then the initial velocity for the given 
mass flow and density was 9.36 cm/s. Finally, we specified the pressure 
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gradient tern as dp/d^ = - 0.278 dynes/an , which is 2.5 times the 
product of the initial density and the initial velocity squared. This 
parameter determines the shape of the channel. 


Several' fluid properties also need to be specified to enable 
the calculation to be carried forward. The set of equations to be 
solved^ Eqs. (5-7) , contain the vescosity y and the ratio of thermal 


conductivity to frozen specific heat, 
to introduce the Prandtl number 




These can be combined 




^C£^=p£ ■ ■ 


Thus we can use y and Pr to express these transport properties. Ex- 
amination of our pure hydrogen calculations of y. Section IV, indicates 
that a reasonable fit for 30 atm, 1000 K < T < 5000 K is 


y = 7.5 E-5 + 1.24 E-7 T (g/cm-s). 


This value was used in the calculation, although it is recognized that, 
as pointed out at the end of Section IV, the viscosity of a mixture 
including cesium will be higher because of its high molecular weight. 
Therefore, the viscosity should be improved in the final version of 
the LHTE program. 

The frozen specific heat can be found from the expressions 
indicated in Section IV for c^, by ignoring the terms involving deriva- 
tives of the number densities. Likewxse can be found by using 

the X^ and X. ^ contributions to X, given in Eqs. (4-55a) and 4-55b) . 
tr int 

When this is done at 30 atm, 1000 K<T < 6000 K, and the values of Pr 
are found, it turns out that Pr varies slowly from 0.68 to 0.67. This 
is in accord with the well-known fact that Prandtl number is a slowly 
varying function of temperature, in the sample calculation being described, 
a constant value of 0.65 was used, though a slightly higher value would 
have been preferable. 
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In order to reduce the complication involved in this first 
calculation, we used Le^ = 1 , thus assuming that diffusion of heat 
and species were the same. It is known that this approximation has 
only a small effect on convective heat flux, compared to using more 
accurate values of Le^, in reacting gasdynamics. 

1 2 

At the initial station, where p and u are constant, ^ ~ "2 ' 

so the value at r is 0.8287, Thus ^ goes from 0 to 0.8287. For this 
w 

first calculation we used only ten intervals in 'f, equally spaced. 

This is not enough to give very accurate results, especially in the 
high radial gradient region near the wall, but we wished only to test 
the calculation scheme and the nature of the results, not produce great 
accuracy. 

The calculation was carried out to ^ = 100 cm, and some of 
the results are presented in the next figures. In Fig. 5-5 are pre- 
sented radial profiles of the laser intensity I at a number of ^ stations.* 
The notable feature is the rapid absorption of laser energy near the 

center of the channel. By 1.5 cm, the center intensity is reduced 

4 2 2 

from 10 to 7000 W/cm , and by 3.5 cm to 1500 W/cm . Meanwhile, there 

is almost no absorption outside a radius of 6 cm. The reason for this 

is the low temperature in the outer regions, where the laser intensity 

is low so the heating is low, and also the wall is constrained to stay 

at 1000 K. 

A similar plot of the temperature profile is presented in Fig. 

5-6. This starts uniform at 1000 K at ^ = 0, but heats rapidly up 
to 6000 K at the center by 1.33 cm. It then begins to cool by radiation, 
and becomes flattened at the center while the hot region spreads outward. 
At 7.1 cm, when Fig. 5-5 shows that most of the laser energy has been 
absorbed, the central temperature is down to *4000 K, and after that 
it continues to decrease. It appears that after the laser energy has 
been mostly absorbed, only energy loss prevails, and one would then 
want to begin the gas acceleration by sharply reducing the channel area. 

* The slope discontinuities in the curves in Figs. 5-5 and 5-6 are caused 
by the plotting routine. The calculated curves have continuous slopes. 
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Fig. 5-5a Laser Intensity Profiles at Various Axial Stations 
0 < ^ < 1.326 C3n. 
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Fig. 5-5b Laser Intensity Profiles at Various Axial Stations 
1.538 < C < 96.529 ati. 
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Actually/ the area is being reduced continuously in this calcu- 
lation, as is clear from Fig. 5-6, where the point at 1000 K is the 
wall at each ^ station. To make this clearer, a plot of ^^.vs ^ is 
given in Fig.'« 5-7. The choice of dp/d^ used produced the channel shape 
shown. By 5 = 10 an the wall radius is down by a factor of 2, and 
by 100 on it is about one quarter of its initial value. This change 
has produced an acceleration, since at 96.5 cm the velocity on the 
centerline is 362 cm/s, compared with its initial value of 9.36 on/s. 

But this is only a small change in velocity, and the real acceleration 
of the flow will need a nozzle to expand it to supersonic speeds, of 
course. Back at 7.1 on, where the laser energy is nearly fully absorbed, 
the centerline velocity is 142 cm/s. Most of the energy has gone into 
heating, not acceleration. 


The absorption of laser energy can best be illustrated by calcu- 
lation of the power remaining in the laser beam at any station C/ by 


r 

/ w 

27Tr I dr . 

0 


i 


This quantity is plotted in Fig. 5-8 against It shows a rapid drop 
beginning at 1 cm. The power is reduced to 10% of its initial value 
by 8 cm, and to 1% by 25 cm. While some of this reduction is due to 
the laser beam being cut off by the shrinking channel size, most is 
due to absorption in the gas. 

The object of the laser heating is, of course, to get the energy 
into the gas. The power is the gas is 

r 

2 

. = / 2Trr (h + u /2) dr 

0 

although for the region calculated the kinetic energy contribution 

is negligible. P is also shown in Fig. 5-8. It rises rapidly as 

G 
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Fig. 5-7 Axial Distr 







the laser energy is absorbed, but losses also rise as the temperature 
increases, so the power in the gas reaches a peak at about 3.6 cm and 
then decreases slowly. The peak power in the gas is e±»ut 0.18 MW, 
which is 18% of the input laser power. At this station the laser power 
is 0.39 MW, so 0.61 MW has been extracted from the beam. If 0.18 MW 
is in the gas, then the remaining 0.43 MW has been lost to the wall, 
either because the outer edges of the beam hit the wall as the wall 
radius decreased, or because of heat flux to the wall. 

The amount cut off by the converging walls can be estimated 

by use of Fig. 5-1. The radius at 3.6 an is 8.1 cm, down from 11.8 

an at the initial station. The P /P, curve of Fig. 5-1 shews that 

r 1 

about 10% of the incoming beam has been cut off by that radius reduction, 
or 0.1 MW, leaving 0.33 MW to be accounted for by wall heat, flux. 

Figure 5-9 shows the distribution of convective and total local 
wall heat flux and where q^ is from Eq. (5-32) and q^ = q^, + q , 
with q^ from Eq. (5-36) . The radiative heating is by far the dominant 

3 

term as can be seen from the factor of 10 difference in the scales. 

2 

The peak heating is about 3.3 kW/cm at 2.4 on. An estimate of the 
integral of q^ up to 3.6 an is 0.25 MW, which is about 0.08 MW short 
of adding up to 1 MW. This 8% error can probably be accounted for 
by the crudeness of the 'I'-grid on which this example was calculated. 

Within that accuracy, we may say the energy balances, either being 
in the gas, in the laser beam, intercepted by the walls, or radiated 
to the walls. 

2 

We may note that the peak heat flux shown in Fig. 5-9, 3.3 kW/on , 
is high, but not unreasonable, and can be handled with fairly conven- 
tional cooling methods. 

This calculation was a preliminary one, meant to show that 
the model produces physically meaningful and reasonable results. It 
shows good absorption at high temperature and very little at low tem- 
perature, at least in a distance of 10 on. It shows that at p = 30 atm 
a large fraction of the laser energy put into the gas is radiated to 
the walls , but even this heat flux can probably be hamdled. 
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5,9 Further Development 

There are a number of further developments that are needed 
to make this bode generally useful. Foremost is the need to generalize 




the radiation model to be applicable to alL useful pressure levels 

and mixture compositions. Second, v/e need to get some experience with 

' [ ( 

riinning the code to learn how to choose dp/d^ so as to maintain a desired . 

^ ^ ' " 

channel shape -whi ch probably should be more"^heirry symm^ifical than 
the present shape. A more refined m^h in ;the 'F coordinate is also ^ 
needed, to fill in more points near the\enter^ With uniform input 
flow and' uniform spacing of points, we start with equal increments 
' in rf ,j; not r. Thus in r the points get closer as r increases. This 
is ^desirable to place more pb/nts in the wall "region. However, with 

Ifi . \ : 

only ten intervals in there is inadequate definition of the central 


region where the large heating is taking, .^place 

I g- - rr - a cju v. - - - - * 




Of ! 


h ^ ll ^ 

Work is proceeding to incorporate these further developments 


into the code to make it generally useful for' parametric studies. 
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VI . CONCLUSIONS 


This interim report describes work which has focused mainly on 
the physical properties of the gases which are important for modeling 
a CW laser-heated thruster. We have described the laser absorption pro- 
perties, the radiation properties, and the thermodynamic and transport 
properties. We have also described a flow code to calculate the channel 
flow of a laser-heated gas, but only one preliminary calculation has 
been made so far. 

The conclusions drawn, then, relate to the gas properties, and 
are based on calculations made with the various gas property models des- 
cribed. 

1. For 10.6 ym radiation, a mixture of 1 part cesium to 99 parts 
atomic hydrogen substantially increases the absorption coefficient at 
low temperatures. The value of 0.1 cm is reached at 2000 K at 100 atm, 
3000 K at 10 atm, and 12000 K at 1 atm. 

2. For 10.6 ym radiation, water vapor has a minimrun absorption 
coefficient at about 500 K, and appears to be quite low below 1000 K. 
However, the values in the vicinity of 500 K are not well-known, and 

the theoretical predictions have not been checked up to 1500 K by experi- 
ments. Experimental measurements between room temperature and 1500 K 
are needed. 

3. For 5.3 ym radiation, water vapor is a good absorber from 
300 to 3000 K. NO is a good low temperature absorber, and will not dis- 
sociate even though it is not an equilibrium state. CO is a good absorber 
above 2000 K, and remains chemically stable up to at least 6000 K. 

4. For 10.6 ym radiation, a mixture of up to 5% water vapor in 
hydrogen does not provide good low temperature absorption. 

5. For 5.3 ym radiation, both NO and water vapor up to 5% in 
hydrogen provide adequate absorption at low temperatures, though NO is 
the better absorber. 
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6. The continuously over-lapping line model of molecular radiation 
provides results in good agreement with other calculations and with experi- 
ments. 

7. Cesiimi is a strong radiator, and its radiation will be an 
important contributor for thrusters of operational size (10 an radius) . 

It is less important for laboratory size thrusters (1 cm radius) , but 
still strong in certain bands. 

8. The thermal conductivity of cesivim-seeded hydrogen is nearly 
the same as for pure hydrogen. The viscosity, however, increases con- 
siderably because of the high molecular weight of cesiian. 

9. The channel flow model of a laser-heated mixture of hydrogen, 
water vapor and cesium produced a preliminary calculation which seems 

to be physically reasonable. 
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